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UCSB South Pole 1994 CMB anisotropy measurement constraints 
on open and flat-A CDM cosmogonies 

Ken Ganga^, Bharat Ratra^, Joshua 0. Gundersen"^, and Naoshi Sugiyama'* 

ABSTRACT 

We develop methods to account for experimental and observational 
uncertainties in likelihood analyses of data from cosmic microwave background 
(CMB) anisotropy experiments and apply them to an analysis of the UCSB 
South Pole 1994 (SP94) experiment. Observationally motivated open and 
spatially-flat A, cold dark matter cosmogonies are considered. Among the 
models we consider, the full SP94 data set is most consistent with Oq ~ 0.1 — 0.2 
open models and less so with old (to ^ 15 — 16 Gyr), high baryon density (0^ ^ 
0.0175/^-2), low density (Oq ~ 0.2 - 0.4), flat-A models. The SP94 data do not 
rule out any of the models we consider at the 2cr level. The SP94 experiment is 
most sensitive to anisotropics on a somewhat larger, model-dependent, angular 
scale than the scale at which the window function peaks. 

For establishing the significance of a detection of CMB anisotropy we derive 
limits using the highest posterior density (HPD) prescription, since it yields 
smaller lower limits. Since HPD limits lead to tighter constraints on the CMB 
amplitude, they also provide for greater discrimination between models. Model 
normalizations deduced from the SP94 data subsets are mostly consistent with 
those deduced from the two-year COBE-DMK data, although the Ka-band 
data prefer a normalization ~ Icr lower than do the Q-band data, the Q and Ka 
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+ Q data favour a slightly higher normalization for the Oq = 0.1 open model 
than does the DMR, and the Ka and Ka + Q data prefer a somewhat lower 
normalization for the older, higher 0^, low-density A models than does the 
DMR. 

Subject headings: cosmic microwave background — cosmology: observations — 
large-scale structure of the universe 

1. Introduction 

Low-density cold dark matter (CDM) cosmogonies, either with flat spatial hypersurfaces 
and a cosmological constant. A, or with open spatial hypersurfaces and no A, are consistent 
with a large fraction of present observational data. For flat-A models, see Peebles (1993), 
Stompor, Gorski, & Banday (1995), Scott, Silk, & White (1995), Ostriker & Steinhardt 
(1995), and Ratra et al. (1997, hereafter RSBG). For open models, see Ratra & Peebles 
(1994), Kamionkowski et al. (1994), Gorski et al. (1995), Liddle et al. (1996), and RSBG. 

These observational data include: 

• measurements of the Hubble parameter Hq (= 100/ikm s~^ Mpc~^), which suggest 

h > 0.55 (e.g.. Pierce & Jacoby 1995; Mould et al. 1995; Tanvir et al. 1995; Hoflich & 
Khokhlov 1996; Baum et al. 1995; Whitmore et al. 1995; but also see, e.g., Nugent 
et al. 1995), and globular cluster age estimates which suggest that the age of the 
universe to > 11 Gyr (e.g., Bolte & Hogan 1995; Jimenez et al. 1996; Chaboyer et al. 
1996); 

• dynamical estimates of the mass clustered on scales ^ 10/i~^Mpc, which suggest 
Oo ~ 0.05 — 0.35 (Peebles 1993) [here the present value of the clustered-mass density 
parameter Oq = 9>TTGph{to) j {?>Ho^), where G is the gravitational constant and pb{to) is 
the mean clustered-mass density now]; 

• estimates of Oq from measurements of the plasma mass fraction of the rich clusters, 
standard nucleosynthesis theory, and the observed light element abundances (e.g.. 
White & Fabian 1995; David, Jones, & Forman 1995); 

• dynamical estimates of the mass clustered on scales ^ 10/i~^Mpc (e.g.. Cole, Fisher, 
& Weinberg 1995; Loveday et al. 1996; Baugh 1996; there is, however, some scatter; 
see, e.g., Dekel et al. 1993; Shaya, Peebles, & TuUy 1995); 
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• the shape of the observed galaxy fluctuation power spectrum (e.g., Peacock & Dodds 
1994; Tadros & Efstathiou 1995); 

• the observed cluster mass and correlation functions (e.g., Liddle et al. 1996; Viana & 
Liddle 1996); 

• the presence of high redshift [z ~ 0.5 — 1) clusters (e.g., Mellier et al. 1994; Luppino 
& Gioia 1995; Small & Dickinson 1995), high redshift [z > 2) galaxy groups (e.g., 
Francis et al. 1996; Pascarelle et al. 1996), high redshift (^^ ^ 3 — 4) damped Lyman-a 
systems (e.g., Giallongo et al. 1994; Lu et al. 1996; Wampler et al. 1996), and the 
observed similarity of the giant elliptical luminosity functions at 2; = and 2; ~ 1 
(e.g., Djorgovski et al. 1995; Glazebrook et al. 1995; Lilly et al. 1995); 

• the large angular scale cosmic microwave background (CMB) anisotropy measured 
by the COBE-DMK experiment and analyzed in the context of these cosmogonical 
models (e.g., Bunn & Sugiyama 1995; Gorski et al. 1995; Stompor et al. 1995); and 

• the smaller-scale CMB spatial anisotropy detected by a number of experiments 
(RSBG; Ganga, Ratra, & Sugiyama 1996, hereafter GRS). 

We emphasize that a number of the above observational estimates are still tentative, 
and that some have large error bars. In particular, any individual set of measurements 
does not have the power to seriously constrain the range of cosmological parameters such 
as Oo and h. It is only the combination of all these different measurements that indicates a 
low-density cosmogony may be a more appropriate model of the universe. 

Following the DMR discovery of large-scale CMB anisotropy and the conflrmation by 
FIRS (Smoot et al. 1992; Bennett et al. 1992; Wright et al. 1992; Ganga et al. 1994), 
attention has been focussed on attempts to use the DMR data to constrain cosmogonies 
which are consistent with other cosmological data (e.g., Bunn & Sugiyama 1995, Gorski 
et al. 1995, Stompor et al. 1995, Yamamoto & Bunn 1996). While the DMR data does 
allow for a much more accurate determination of the normalization of cosmological models 
than has previously been possible, it does not have the discriminative power to signiflcantly 
constrain other cosmological parameters. 

Besides the DMR and FIRS measurements of large-scale CMB anisotropy, there are 
now a number of other detections of anisotropy on smaller angular scales (Hancock et al. 
1996; Piccirillo et al. 1997; Netterfleld et al. 1997; Gundersen et al. 1995, hereafter G95; 
Ruhl et al. 1995; de Bernardis et al. 1994; Tanaka et al. 1996; Cheng et al. 1996; Griffin 
et al. 1997). It is of some interest to determine whether these smaller-scale anisotropy 
measurements can help to constrain cosmological model parameters. Given the error bars 
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associated with these measurements, it is not surprising that, even if one makes use of the 
DMR data to fix one free parameter and normahze cosmological models, all the smaller-scale 
anisotropy experiments combined do not very significantly constrain other cosmological 
parameters (RSBG; GRS). In particular, a quantitative goodness-of-fit comparison of the 
CMB anisotropy predictions from two-year DMR-normalized, open and fiat-A CDM models 
to all available CMB detections by GRS found that in most cases the DMR error bars 
precluded robust conclusions about model viability.^ 

The comparisons between the theoretical predictions and the data done in RSBG and 
GRS did not make use of the complete data from each experiment. Rather, they used a 
single number with error bars related to the rms anisotropy detected by the experiment. 
There are two significant issues related to such a sparse representation of the observational 
data from an experiment. Since the procedure used to extract this single number from 
the observational data set typically assumes a fiat bandpower spectrum or a gaussian 
autocorrelation function for the anisotropy, it is of interest to examine how sensitively the 
final result depends on this assumption and to determine the result using more realistic 
forms for the spectrum. Secondly, although the experiments are sensitive to a range of 
spatial scales, representing the observational data by a single number effectively discards 
additional spatial correlation information. It is therefore also of interest to examine whether 
this additional spatial correlation information aids in discriminating between model spectra. 

In this paper we address these issues for the SP94 data set of G95, using anisotropy 
spectra derived in gaussian, adiabatic, open and flat-A, CDM cosmogonies (RSBG). Bond 
& Jaffe (1996) have also analyzed the SP94 data, mostly making use of CMB anisotropy 
spectra in spatially-flat models, but considering a wider variety of these models than we 
have here. The SP94 windows are sensitive to anisotropy on angular scales where the 
open and flat-A model CMB spatial anisotropy spectra are not very scale dependent, so 
we expect only a weak dependence of the deduced signal amplitude on the assumed CMB 
anisotropy shape. 

SP94 is the most recent of the SP CMB anisotropy experiments at the South Pole. 
The observational results of SP89 and the ACME telescope are discussed by Meinhold & 



®The new four-year DMR results (e.g., Bennett et al. 1996; Gorski et al. 1996) indicate a normalization 
slightly lower than that found from the two-year data. This is mostly a consequence of more detailed 
modelling of foreground Galactic emission (Banday et al. 1997; Kogut et al. 1996). From Figs. 5 and 6 
of GRS, which correspond to the models normalized at the —la value of the two-year normalization (or, 
^ cr/3 below the nominal value of the four- year DMR normalization), one sees that most models considered 
by GRS and here are consistent, in the goodness-of-fit sense, with the small-scale CMB anisotropy detection 
data. 
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Lubin (1991) and Meinhold et al. (1993). Analyses of the SP89 data, in the context of 
fiducial scale-invariant, spatially-fiat CDM models and some isocurvature models is given in 
Vittorio et al. (1991) and Bond et al. (1991). The SP91 observational results are in Gaier 
et al. (1992) and Schuster et al. (1993). The Gaier et al. data have been analyzed in the 
context of fiducial CDM by Gorski, Stompor, & Juszkiewicz (1993), Muciaccia et al. (1993) 
(who also consider tilted CDM models), Dodelson & Jubas (1993), Bunn et al. (1994), 
and Stompor & Gorski (1994). Stompor & Gorski also consider the now more fashionable 
fiat-A CDM models and analyze the Schuster et al. (1993) data. These analyses only used 
the data from the highest-frequency channel of the SP91 scans to draw conclusions about 
the CMB anisotropy. Off-diagonal channel-channel noise correlations were accounted for in 
the multifrequency reanalysis of SP91 given in G95, who made use of the flat bandpower 
approximation. The inclusion of these off-diagonal correlations has little effect on the most 
probable sky signal amplitude. However, they increase the iblcr error bars by ~ 20%. 

The SP94 experiment extends the SP91 frequency range by taking data with a Q-band 
receiver, as well as with a Ka-band receiver. Descriptions of the SP94 experiment may be 
found in Gundersen et al. (1994), G95, and Gundersen (1995). The far-fleld beam patterns 
of the Ka-band and Q-band radiometers, coupled to the ACME telescope, were measured in 
both azimuth and elevation, at 27.7 GHz for Ka and 41.5 GHz for Q. These measurements 
indicate that in both azimuth and elevation the beams are well-approximated by gaussians 
down to 30 dB (Gundersen 1995). The beamwidths measured in azimuth and elevation 
differ from the average beamwidth by ~ 1.5 — 2% (Icr). In our analysis, we do not explicitly 
account for this ellipticity but assume circular beams with beamwidth errors bars that do 
account for this difference. Multifrequency beam patterns estimated from observations of 
the Moon with the Q-band system and observations of the Eta Carina region with both 
systems are consistent with the far-fleld beam measurements. Ear-fleld multifrequency 
two-dimensional beam maps of the far side lobes have not been made. Based on the 
one-dimensional far-fleld measurements of the beam patterns, the main beam efflciency is 
estimated to be > 98%. The < 2% loss to the beam side lobes is not corrected for in this 
analysis. A more detailed discussion of these issues may be found in Gundersen (1995). 

The Ka-band (26 — 36 GHz) is multiplexed into four channels centered at u = 27.25, 
29.75, 32.25, and 34.75 GHz, with 3 dB bandwidths of 2.5 GHz. The measurements 
described above indicate that the individual Ka channels have a frequency-dependent 
gaussian beamwidth aff^^ = (0.70 ±0.04)° x [27.7 GHz/i/] (1 cr error). The Q-band (38 — 45 
GHz) is multiplexed into three equal-width channels centered at u = 39.15, 41.45, and 43.75 
GHz with 3 dB bandwidths of 2.3 GHz and a frequency-dependent gaussian beamwidth 
"^l?^ = (0.47 lb 0.04)° X [41.5 GHz/i/] (Icr error). The measured passband central frequency 
errors are ~ 1%, and are ignored in our analysis. Ear-fleld measurements of the Ka- and 
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Q-band beamwidths were performed at 27.7 and 41.5 GHz respectively. We therefore 
assume that all four Ka-band beamwidth uncertainties and all three Q-band beamwidth 
uncertainties shift together. However, this is not necessarily the case, since the Ka and 
Q measurements were performed consecutively, with different horns, but with the same 
telescope optics. 

The SP94 Q and Ka observations were performed consecutively. Since the multiplexing 
was done after amplification and initial filtering, the HEMT amplifiers and the atmosphere 
introduce significant intraband channel-channel correlations, which are accounted for in the 
analysis. 

Unlike the stepped scans of SP89 and SP91, SP94 data were taken during smooth, 
azimuthal, constant declination, constant velocity scans extending 20° on the sky. While 
observing, the beam was sinusoidally chopped with a half peak-to-peak chop amplitude 
of 1.5° on the sky. The rms fractional uncertainty in the peak-to-peak chop angle was 
measured to be 0.04%, and is ignored in our analysis. The two-beam/single-difference data 
were binned into 43 bins, with a bin size of (20/43)° on the sky. As discussed in G95, the 
relative-pointing uncertainty "on the sky" is ±0.02° in azimuth and ±0.05° in elevation, 
and is ignored in our analysis here. The absolute-pointing uncertainty on the sky is ±0.12° 
in azimuth and in elevation (G95). The absolute-pointing uncertainty is irrelevant for either 
the Ka- or Q-scan data analyses, but could be a significant issue when the Ka and Q data 
are combined for the Ka ± Q analysis, or when comparing the results of the Ka and Q 
analyses. However, for reasons discussed below, this is ignored in the analysis here. 

The G95 data were taken in smooth scans of ±10° on the sky centered at a = 45° and 
8 = —61.8° (1994), with a total sky coverage ~ 20 deg^. The data editing, as described in 
G95, resulted in the removal of ~ 25% of the data — less for Ka, more for Q — and left 88 
h and 108 h of usable Ka- and Q-band data respectively. After editing, an offset and linear 
gradient were removed from the data for each channel and each scan (of angular extent 
20° on the sky; we note that the SP group sometimes calls this a half-scan), and for each 
channel the binned data from the individual scans were coadded (G95). The data were also 
corrected for atmospheric absorption (G95). This is the reduced data we analyze in this 
paper. 

While in this analysis we assume that the SP94 data is purely CMB anisotropy, it 
is prudent to bear in mind that one cannot yet conclusively rule out a small amount of 
undetected, non-CMB contamination in the data. The G95 time subset check, which 
shows that when the data is divided into 4 contiguous time series, all 4 subsets and 
all combinations of 3 of the 4 subsets are mutually consistent, indicates that transient 
atmospheric phenomena are not a significant issue. G95 also show that the contamination 



- 7 - 



expected from Galactic diffuse synchrotron emission, from the Sunyaev-Zel'dovich effect, 
and from Galactic 20 K dust emission are small compared to the signal they see. See 
Gundersen (f995) for a more detailed discussion. G95 could not, however, rule out 
non-CMB foreground discrete radio source contamination (also see Gundersen et al. f996). 

The Ka-band and Q-band radiometers were calibrated using a combination of different 
methods. Without the ACME telescope, the radiometer V/K were determined for an 
ambient temperature Eccosorb load and for a liquid nitrogen cold load. These were used in 
combination with sky zenith scans and an atmosphere model to determine the atmospheric 
contribution to the antenna temperature and to derive a sky zenith temperature. The 
sky zenith scans were performed with the radiometers coupled to the ACME telescope. 
The calibration of the radiometers, determined from the three different combinations of 
ambient load, cold load, and sky load, were found to have a channel to channel uncertainty 
of 3% and an absolute uncertainty of fO% (Ic). This determination of the calibration was 
consistent with that determined from a cryogenic termination that was used to produce 
a load similar to that of the sum of the atmosphere and the CMB. Measured microwave 
emission from the Moon was also consistent with a model of the Moon's emission (Keihm 
f983) to within the ~ 20% uncertainty. Some CMB anisotropy experiments calibrate on a 
point source to determine V/Jy. SP94 does not, so for SP94 the beamwidth uncertainty 
does not influence the calibration uncertainty. In our analysis, the SP94 beamwidth and 
radiometer-calibration uncertainties are treated as independent sources of uncertainty. 
A more detailed discussion of these points may be found in Gundersen (1995). In what 
follows, we treat the absolute-radiometer-calibration uncertainty as an additional, purely 
statistical, uncertainty to be included in the likelihood analyses. 

Besides accounting for beamwidth and calibration uncertainties, and using anisotropy 
spectra in observationally motivated cosmogonies, our analysis differs from that of G95 in 
two other respects. Unlike G95, we do not make use of data-weighted windows to quote 
limits, and we also use a different statistical prescription to determine limits from the 
probability density distribution functions. A more detailed discussion of these points is 
given in §2 below. 

In §2 we outline the computational techniques used in our analysis. Our results and a 
discussion are in §3, and we conclude in §4. 
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2. Computation 

The CMB spatial anisotropy temperature variations, 5T/T^ can be decomposed into 
spherical harmonics, 

7f^ = 2^ 2^ aimYl.ra{0,(p). (1) 

1=2 m = -l 

For gaussian models we may characterize the anisotropy by its power spectrum, C;, which 
is defined in terms of the ensemble average 

(«/m«*'m') = Cl^ll'^mm'- (2) 



One would like to use anisotropy observations to measure the power spectrum over the 
entire range of / to which the experiment is sensitive. However, the broad width of the SP94 
windows and the limited sky coverage prevents one from extracting each individual C;. It 
is therefore necessary to assume a functional form for the C;s over the range of / which 
SP94 is sensitive to, with the overall normalization allowed to be a free parameter and with 
the shape of the C; allowed to depend on a small number of other parameters. One may 
then compare these C;s to the spectrum of the SP94 data and determine the value of the 
normalization and of the shape parameters which best reproduce the data spectrum. 

On large angular scales, when the effects of the pressure of the photon-baryon fluid 
and matter velocity perturbations at photon decoupling can be ignored, an approximation 
to the CMB anisotropy spectrum in the flducial CDM model is provided by the flat CMB 
angular spectrum (Peebles 1982), 

^ 6C2 ^ 24^ (Qrms-Ps/ro)' . . 

' + 5 + ' ^' 

where Qrms-PS is the corresponding quadrupole-moment amplitude of the model CMB 
anisotropy and To is the CMB temperature now. This spectrum, normalized to best 
reproduce the two-year DMR sky maps is shown in Figure 1 (line labelled "Flat"; Gorski et 
al. 1994). The DMR range for Qrms-PS in this model is given in the last line of the third 
column of Table 12. 

The flducial CDM model is an Finstein-de Sitter model with gaussian, adiabatic, 
scale-invariant energy density perturbations (Harrison 1970; Peebles & Yu 1970; Zel'dovich 
1972), standard recombination (Peebles 1993, §6), Hubble parameter h = 0.5, and baryon 
density parameter 0^ = 0.0125/i~^. Scale-invariant energy density perturbations are 
generated by quantum-mechanical fluctuations during an early epoch of inflation in a 
variety of spatially-flat inflation models (e.g., Fischler, Ratra, & Susskind 1985). In inflation 
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models the small observed CMB anisotropy could be the consequence of the small ratio of 
the inflation epoch mass scale to the Planck mass (Ratra f99f , and references therein; also 
see Banks et al. f995). 

The spectrum for this model (Of4 in Table 2), normalized to best reproduce the 
two-year DMR sky maps, is shown in Figure f (Gorski et al. f995). The range for Qrms-PS 
allowed by the DMR and a subset of the SP94 data is given in the line labelled 014 of 
Table 12. Note that Models 014 and A12 have identical CMB anisotropy spectral shape, 
but are normalized using the slightly different DMR galactic- and ecliptic-coordinate maps. 
The fiducial CDM model is no longer thought to provide an adequate representation of the 
observed universe. 

On the other hand, low-density open and fiat-A CDM cosmogonies are consistent 
with present observational data. For historical reasons (Guth 1981, also see Kazanas 
1980; Sato 1981a,b), no longer vahd (Gott 1982; Guth & Weinberg 1983), the low-density 
A CDM models considered here are taken to have fiat spatial sections, so the simplest 
power spectrum for gaussian, adiabatic, energy-density perturbations in these models is 
the scale-invariant one. Such a spectrum is generated by quantum zero-point fiuctuations 
during an early epoch of infiation in a spatially-fiat infiation cosmogony. The simplest 
spectrum for gaussian adiabatic energy-density perturbations consistent with open spatial 
sections is that generated by zero-point fiuctuations during an early epoch of infiation in 
an open model (Ratra & Peebles 1994, 1995; Bucher, Goldhaber, & Turok 1995; Lyth & 
Woszczyna 1995; Yamamoto, Sasaki, & Tanaka 1995; Bucher & Turok 1995). 

The values of the parameters Oq, ^, and 0^ used here are chosen to be roughly 
consistent with present observational estimates of Oq, ^, the age of the universe, and 
the constraints on 0^ that follow from the observed light element abundances in the 
standard nucleosynthesis model (RSBG). Two standard-recombination open and two 
standard-recombination fiat-A spectra are shown in Figure 1. 

In this analysis we wish to draw conclusions about parameters representing effects that 
must exist; i.e., Oq, ^, and 0^. To simplify the comparison to SP94 data, the effects of 
tilt, primordial gravity waves, and reionization are ignored. Tilt, gravity waves, and early 
reionization are unlikely to be significant in viable open models, although there are mild 
indications that some such effect might be required to bring some flat-A CDM models into 
agreement with observational data (Stompor et al. 1995; Scott et al. 1995; Ostriker & 
Steinhardt 1995; RSBG; Klypin, Primack, & Holtzman 1996; GRS). We emphasize that 
any effect which significantly modifies the shape of the spectra used in our analysis here, 
such as isocurvature perturbations or nongaussian anisotropy, will affect the final results. 
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The general procedure for the computation of the CMB spectra is discussed by 
Sugiyama (1995). A discussion of the various assumptions and approximations involved in 
such a computation, as well as a discussion of the resulting accuracy of the computation, 
may be found in Hu et al. (1995). 

Figure 2 shows the zero-lag windows, W;, for the seven individual SP94 channels, at 
the nominal beamwidths. Figure 2 also shows the uncertainties due to the Icr beamwidth 
uncertainty in the highest and lowest frequency channels. 

It has become conventional to summarize anisotropy experiment windows in terms of a 
few parameters (e.g.. Bond 1996). These are the value of / where Wi is greatest, /m, the two 
values of /, /e-^ s, where Wi = e~'^'^Wi^^ and the effective multipole /g = I (IWi) / 1 (Wi) ^ 
where 

The values of these parameters for each of the seven channels are given in Table 1, both 
for the nominal beamwidths, and for beamwidths Icr smaller and larger than the nominal 
ones.^ 

The range of multipole moments to which an anisotropy experiment is sensitive 
depends on both the window Wi and the sky signal C;. Depending on the form of the sky 
signal, the window parameters /g, ^m, and /g-0.5 do not necessarily give a good indication of 
the multipoles to which an experiment is sensitive. 

CMB anisotropy experiments are sensitive to ^Trms, the rms of the anisotropy seen 
through their window. Here 

00 

{ST^msf = E(<^r™s')/, (5) 
1=2 

where 

[ST^ms')! = To'^-^^CWi, (6) 

and To = 2.726 it 0.010 K (Mather et al. 1994). Given an anisotropy model, (^Trms^)/ 
provides a convenient way of establishing the approximate range of multipoles / an 
experiment is sensitive to. It is approximate since it ignores spatial correlation information. 
In Figure 3 we plot (^Trms^)/ for the models shown in Figure 1, for two SP94 window 
functions. In Table 2 we list /m, the value of / at which (^Trms^)/ is at a maximum, and /g-0.5 



®Note that in our internal computations these and other numerical values are not truncated at 2 or 3 
significant figures as in the tables; as a consequence, when the truncated numerical values of the tables are 
used to rederive some of our results there will be small differences from the values listed in the tables. 
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the two multipoles where (^Trms^)/ _o 5 = ^ '^'^{ST,-ms'^)i^^ for these two windows, and for all 
the models we consider here. 

The range of multipole moments to which each channel is sensitive is quite model 
dependent, and this range is not adequately summarized by the Wi parameters /g and 
/e-0.5 . For the spectra considered here, the SP94 experiment is sensitive to a somewhat 
larger angular scale than the effective angular scale determined by the Wi parameter 
/g. In this case, however, the (^Trms^)/ parameters 1^ and /g-0.5 do provide an adequate, 
model-dependent, characterization of the range of / to which the SP94 experiment is 
sensitive. This is not true in general, and only works here because the models considered 
here have spectra that are fairly smooth in the relevant range of /-space. 

In what follows, we shall have need for the bandtemperature (e.g.. Bond 1996), 

ST, = 4^^. (7) 



Note that despite the notation, 5Ti does not depend explicitly on /. Compared to ^Trms, ^Ti 
has the advantage of being insensitive to the normalization of Wi. As we shall see below, it 
is still fairly sensitive to the precise shape of Wi. 

The reduced SP94 data used in the analysis of G95 are shown in Figure 4. In the 
line labelled "Sky" in Table 3, we give the estimated rms of the anisotropy, computed 
from the data of Figure 4 as the square root of the difference between the variance of the 
mean temperatures and the variance of the error bars. As this ignores off-diagonal noise 
and spatial correlations, as well as the oversampling of points on the sky'', which are all 
accounted for in the likelihood analysis, it could deviate from the true sky rms. 

For the purpose of the following discussion we assume that the SP94 sky signal is 
purely CMB anisotropy without any offset or gradient, which would have been removed 
while fitting out an experimental drift. Following Bond et al. (1991), the likelihood function 
at a given beamwidth and nominal calibration (hereafter, the "bare" likelihood function) 
for a given model C; is computed according to 

Lid) oc ^=^^^^_e->^^/2, (8) 
/det(FM-if) det(M) 



where 



X 



A^(M-i - M-if (FM-if)-iFM-i )A. (9) 



"^The SP94 FWHM beamwidths are significantly larger than the data-bin-separation on the sky. 
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Here A is the vector of temperature deviations (with = 43 X 7 elements for the full data 
set, for example) and M = C + S is the correlation matrix of the data in the context of the 
model considered. C represents correlations arising from the specific model considered and 
can be computed using eq. (2) of G95, while S represents instrumental and atmospheric 
noise. The unwieldy term in the expression for arises from assuming a uniform prior in 
the amplitude of the drifts and marginalizing over all possible values of offset and gradient, 
f is the array of functions that have been fit out of the data, in this case an offset and 
gradient from each channel; it is thus an X 2 matrix for the full data set. Here, because 
the bin weights for each channel are fairly uniform, the difference between using unweighted 
fitting functions and weighted ones, as recommended by Bunn et al. (1994), is small. This 
would not be the case, however, for experiments such as FIRS which do not uniformly 
sample the area covered. As discussed above, we parametrize the model C; using Qrms-PS, 
Oo, and 0^. 

With this prescription, for a given data set, all models yield the same likelihood at 
Qrms-PS = /iK, i.e., no sky signal. This provides a convenient way to normalize across 



In Table 3, the lines labelled "FBP" are central ^Trms values derived from likelihood 
analyses using the fiat bandpower (FBP) spectrum (eq. [3]) at the nominal value of the 
beamwidths for each channel. 

To derive the Qrms-PS central value and limits from the likelihood function for a given 
Ci and data set, we adopt the following prescription (Berger 1985, p. 140; Myers, Readhead, 
& Lawrence 1993, §3.1; Gorski et al. 1996). 

Bayes's theorem states that the posterior probability density distribution is proportional 
to the likelihood function multiplied by the prior probability. Lacking good information 
to the contrary, here we assume a uniform prior in (5rms-ps(^ 0), resulting in a posterior 
probability density distribution equal to the likelihood function. 

The central value is taken to be the value of Qrms-PS at which the probability density 
distribution peaks. The 2(7 limits, Q-2a and Q+2a^ are defined as the two values of Qrms-PS 



and for which (5+2(7 — Q-2a is minimized. This is the highest posterior density prescription, 
hereafter HPD. If Q-2a is greater than /iK, we say that there is a (2(t) detection; we then 
determine Q+ia and Q-ia analogously. If instead Q-2a is zero, we say that there is no 
detection and integrate the probability density function over Qrms-PS, starting from /iK 
until we get to the value of Qrms-PS that includes 1 — 0.5(1 — 0.9545) = 97.72% of the total 



models. 



for which 




(10) 
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area under the probability density function and call this value of Qrms-PS the 2cr upper 
limit. This is the equal tail prescription, hereafter ET. Of course, the choice of how to 
define limits depends on what one wishes to use the data for, and will hopefully not be a 
significant issue when the data improves. The definition we choose to adopt yields smaller 
lower limits for detections and yields larger upper limits for nondetections. It also leads to 
the tightest constraint on Qrms-PS and therefore provides for greater discrimination between 
models. 

In all cases we computed limits up to ±3(7 (99.73% HPD). We do not quote them here, 
but emphasize that even accounting for calibration and beamwidth uncertainties, the Ka, 
Q, and Ka + Q data subsets show evidence for at least 3(7 detections. 

G95 quote ±1(7 ET limits for the Ka, Q, and Ka ± Q data. Their prescription differs 
from that adopted here. In Table 4 we give the G95 central values and ±1(7 ET limits 
along with our central values and ±1(7 limits computed using both the G95 ET prescription 
and the HPD prescription. These are given in terms of bandtemperature for the combined 
Ka, Q, and Ka ± Q data, extracted using the flat bandpower angular spectrum. They 
are calculated ignoring calibration and beamwidth uncertainties. The last two columns of 
Table 4 give the average of the ±1(7 error bars in /iK and as a percentage of the central 
value. 

We now want to account for the additional uncertainties induced by the radiometer- 
calibration and beamwidth uncertainties. 

Most experiments quote a Icr uncertainty in the widths of their beams. We assume 
that the probability distribution for the beamwidth, 6, is gaussian, 

P{b) = ^^e-(^-^o)V(2..^)^ (11) 

V 27r (76 

for b > 0, where bo is the nominal beamwidth and erf, is the uncertainty in the beamwidth. 
The likelihood function marginalized over beamwidth uncertainty is then 

/•oo 

iv(grms-Ps) = / dbP{b)L{Q,^,_ps,b), (12) 
JO 

where iv((5rms-PS, is the bare likelihood function derived if a beamwidth of b is assumed. 
To evaluate the expression in eq. (12), we determine a select few beamwidths at which 
to compute the likelihood function and use Gauss-Hermite quadrature summation to 
approximate the integral. In this case we used three point quadrature (e.g.. Press et al. 
1992), which will be accurate so long as erf, <^ bo. 

As discussed above, we assume that the beamwidth uncertainties for all channels shift 
together. This is justifled for the individual Ka channels and for the individual Q channels. 
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However, since the Ka-band and Q-band measurements were performed consecutively, 
using different horns and the same telescope, the Ka and Q beamwidth uncertainties do 
not necessarily shift together. Since we have analyzed this data set assuming that the Ka- 
and Q-band beamwidth uncertainties shift together, in this case our limits are somewhat 
overconservative. Our neglect of the absolute-pointing uncertainty when combining the 
Ka- and Q-scan data sets (for the Ka -|- Q analysis) compensates approximately for this 
effect. We note that a proper accounting of the absolute-pointing uncertainty would require 
analyses of the data with the Ka- and Q-scan bins offset on the sky. These need to be 
performed at various values of the offset, and one then "marginalizes" over the offset. Other 
uncertainties, such as that in the chop amplitude, which is exceedingly small for SP94 and is 
ignored here, may also be accounted for by a straightforward generalization of our method. 

In passing, we note that the technique used by G95 to measure their beamwidths 
simplifies the issue of accounting for the beamwidth uncertainty. If the beamwidths of the 
different channels are measured at different frequencies, or were found not to factorize into 
a frequency-dependent term and a frequency-independent term, it would be a slightly more 
complex matter to properly account for the beamwidth uncertainty. 

Most experiments quote a Icr fractional uncertainty in their calibration. We assume 
that the probability distribution for the calibration, C, is gaussian: 



for C > 0, where ac is the Icr fractional uncertainty in the calibration, 0.1 in this case, and 
C is greater than 0. 

If C is not equal to 1, then the "real" value of Qrms-PS is related to the "bare" value of 
Qrms-PS, which wc wiU hereafter refer to as Qrms-PS? through 
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g-(C-l)V(2ac^)^ 



(13) 
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The likelihood function marginalized over calibration uncertainty is then 




(15) 



With eqs. (13) and (14), this implies 
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In passing we note that the factor of Qrms-PS in the denominator of the prefactor 
on the right-hand side of this equation makes the hkehhood marginahzed over Qrms-PS 
logarithmically divergent at large Qrms-PS- Since this occurs at extremely large values 
of Qrms-PS, this is uot an issue that need bother us here, and is most likely due to the 
invalidity of our assumption (in eq. [13]), that there is a nonzero probability for C = 0. 

Some intuition about eq. (16) may be gotten by studying how it transforms simple 
analytical "bare" likelihoods. If there is no noise and sample variance, and the "bare" 
likelihood indicates a perfect, Dirac delta, detection at Qo, 

-^(Qrms-Ps) = <^(Qrms-PS ~ Qo 

then from eq. (16), 

HQ rms — PS ) 

i.e., the calibration uncertainty skews, and broadens the "bare" likelihood, and iv((5rms-ps) 
has an amplitude (Qrms-ps) dependent width. If on the other hand, one has a good "bare" 
detection at Qo, for which we assume a gaussian form for the "bare" likelihood, 

\/ Ztt (To 

then in the limit when ctq <^ Qq (and assuming that ctq ccQrms-PS and that 
ccQrms-PS ^ Qo) it may be shown that eq. (16) results in 

iv(grms-Ps) = ^ e-(Qrn.s-PS-Qo)V[2{(^cQ.n.s-Ps)H^on]. (20) 

^27r{(crc(5rms-Ps)^ + C^O^} 

i.e., the calibration uncertainty again skews and broadens the "bare" likelihood function, 
and the calibration-uncertainty-corrected likelihood has an amplitude-dependent width. 
This is what is intuitively expected, and differs a little from the usual naive prescription for 
accounting for the fractional calibration uncertainty by adding it in quadrature to the error 
bars derived from the "bare" likelihood function (e.g., RSBG). 

Although the beamwidth and calibration probability distributions are not known for 
any experiment, the gaussian assumptions (eqs. [11] & [13]) are probably quite reasonable, 
and we use these in what follows. It is important, however, to bear in mind that this is an 
assumption which could be checked by direct measurements. 

Finally, we estimate the anisotropy spectral index (3 (WoUack et al. 1993) using both 
the estimated SP94 "sky" rms, and a flat bandpower model spectrum (eq. [3]) with the 



-(Qr 



-Qo)V[2(<TcQr 
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SP94 Wi. For a frequency spectrum STrmsii^) oc , we have 




(21) 



For a flat bandpower model spectrum, this may be rewritten as 



(3 = In 




(22) 



where I{Wi) is defined in eq. (4), and we have used eq. (7) and the fact that for a fiat 



bandpower angular spectrum STi = a/12/5 Qrms-PS- Such a computation makes sense 



for SP94 since the bandpass filters are fairly well modelled as narrow tophats. The 
central values of (3 are given in Table 5 for the nominal beamwidth SP94 windows (these 
computations use the values for u and I{Wi) for the nominal Wis given in Table I, and 
the values in the row labelled "Sky" of Table 3). The (3 values estimated using eq. (22) 
do not account for off-diagonal noise and spatial correlations, nor for the oversampling of 
points. 

In Table 6 we give bandtemperature central values and both FT and HPD iblcr limits 
derived from flat bandpower likelihood analyses of the individual-channel SP94 data sets. 
These ignore beamwidth and calibration uncertainty. The last two columns of Table 6 give 
the average of the iblcr error bars, in /iK and as a percentage of the central values. 

In Table 7 we give the corresponding numbers derived from flat bandpower likelihood 
analyses of the combined Ka, Q, and (full) Ka -|- Q data subsets, now accounting for 
calibration and beamwidth uncertainties using two prescriptions: (I) that suggested by 
one of us (JOG) to account for these uncertainties by adding 15% of the central value in 
quadrature to the G95 FT error bars, which is the de facto "standard" prescription used 
in the CMB anisotropy held; and (2) that described above to flrst correct the likelihood 
functions for these uncertainties, and then to use these corrected likelihood functions to 
derive the relevant numbers. 

Central values and limits for Qrms-PS, both ignoring and accounting for beamwidth 
and calibration uncertainties, are given in Tables 8 and 9 for the Ka2, Q2, Ka, and Ka -|- Q 
data sets. Some of the likelihood functions used in the derivation of these numerical values 
are shown in Figs. 5. In Fig. 6 we compare the calibration- and beamwidth-uncertainty 
corrected likelihood functions for the selected models. These data sets were chosen for 
illustrative purposes only; the Ka2 data does not have a detection, and the Q2 data has one 
of the best detections. 




In Tables 10 and II we list the values of the probability density distribution functions 
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at the peak, and the marginahzed (over Qrms-ps) probabihty density distribution values®, 
both ignoring and accounting for beamwidth and calibration uncertainties, for the Ka + 
Q, Ka, Ka2, and Q2 data sets. The models shown in the Figures were chosen on the basis 
of their marginal probability distribution values for the Ka + Q data set, accounting for 
beamwidth and calibration uncertainties. For all the combined SP94 data sets, and some of 
the individual-channel ones as well, the likelihood functions are peaked, and well-separated 
from /iK. See Figs. 5 and 6. For instance, for the Ka -|- Q data set the amplitude of the 
likelihood function at the peak is approximately 5 • 10^^ times larger than that at /iK. For 
such likelihood functions it makes sense to choose between models on the basis of the value 
of the marginal probability distribution function. Model 01 (an Oq = 0.1 open model) is 
the most likely model, followed, amongst our selected models, by Flat (flat bandpower), 
models Oil and A2 (an open Oq = 0.5 model and a spatially-flat Oq = 0.2 model), 014 
(fiducial CDM), and, the least likely one, AlO (a fiat-A Oq = 0.4 model). These selected 
models include the most and least likely open and fiat-A ones amongst the models we 
consider here. 

Since both the open and fiat-A model CMB spatial anisotropy spectra shapes depend 
on Oo, and 0^, the 14 open and 12 fiat-A model spectra we use in the analyses here do 
not provide enough resolution in model-parameter space for us to be able to marginalize 
over these parameters to extract most likely values for each separately. 

In Tables 12 and 13 we list central values and limits for Qrms-PS for all the SP94 data 
sets. These account for beamwidth and calibration uncertainties. These tables also give the 
corresponding two-year DMR results, accounting for both statistical and systematic DMR 
uncertainties (Stompor et al. 1995, also see RSBG). Table 14 shows ^T;, computed from 
the Qrms-PS values of Tables 12 and 13 using eqs. (5) - (7). Fxcept for the fiat bandpower 
model, it is not possible to compute 5Ti for the Ka, Q, and Ka -|- Q data sets (since they 
are combinations of data from different windows). In Tables 15 and 16 we list the values 
of the probability density distribution functions at the peak and the marginal probability 
density distribution function values, for all the data sets, accounting for beamwidth- and 
calibration-uncertainty corrections. 

Figures 7 show selected model spectra normalized to the SP94 (Ka, Q, and Ka -|- Q) 
and DMR data sets, as well as the corresponding individual-channel SP94 results. 



®The marginalized probability density distribution function is the probability density distribution function 
for just the cosmological parameters SIq, h, and Sl^. In what follows we usually call this the "marginal 
probability distribution function" . 
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3. Results and Discussion 

From Table 3 we see that the values of the rms estimated from likelihood analyses using 
flat bandpower spectra agree quite well with the value of the "sky" rms estimated from the 
data. For four channels they agree to better than 5% and for six channels they agree to 
better than 25%. This gives some indication of the importance of the off-diagonal noise, 
spatial correlations and oversampling of points, and provides semiquantitative confirmation 
of the likelihood method results. The Ka4 channel "sky" rms is ~ 50% larger than that 
derived from the likelihood analysis. The data for this channel does look somewhat different 
from that for the other channels [see the lowest panel of Fig. 4(a)]. To examine whether 
this can significantly affect the results deduced from the likelihood analyses, we have 
reanalyzed the combined Ka and Ka + Q data subsets with the Ka4 channel excluded from 
the analyses. When the Ka4 channel data is dropped from the Ka analysis, the central 
Qrms-PS values rise by ~ 3 — 6% (depending on model), and the Icr error range broadens 
by ~ 5 — 8%, depending on the model. When it is excluded from the Ka + Q analysis, the 
central values drop a little and the Icr error ranges slightly tighten. These changes are not 
big enough to warrant dropping the Ka4 channel data, so the numbers we quote here are 
from the full analyses, including the Ka4 data. 

In Table 4 we compare the central values and iblcr limits on bandtemperature derived 
from likelihood analyses with the flat bandpower spectrum. The numbers are the 16.0% 
and 84.0% FT hmits derived by G95, and both the 15.87% and 84.13% FT as well as the 
68.27% HPD limits found here for the Ka, Q, and Ka + Q data sets ignoring beamwidth 
and calibration uncertainty. The computations described here use more accurate window 
functions than those used by G95; they differ by ~ 0.5 — 2%. Scaling from analyses done 
at three different beamwidths, we estimate that these window function differences and the 
slight differences in FT limit deflnitions are probably responsible for most of the difference 
between the 5Ti values computed here and in G95. 

Table 4 also gives the 68.27% HPD itlcr limits. We adopt this prescription to quote 
limits for detections, since it results in a more conservative —Icr limit than does the 
corresponding FT prescription. It also leads to slightly tighter average constraints on the 
bandtemperature, and so provides for slightly greater discrimination between models. Both 
these effects are evident in the numbers of Table 4: comparing, e.g., the third and fourth 
rows, we see that the HPD prescription results in somewhat smaller upper error bars and 
somewhat larger lower error bars than does the FT prescription. 

Unlike the Saskatoon (SK) experiments (WoUack et al. 1993; Netterfleld et al. 1995, 
hereafter N95; Netterfleld 1995; Netterfleld et al. 1997) where the shape of the window 
functions are frequency-independent, the SP94 window functions are frequency-dependent. 
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As this couples spatial (/) and frequency (u) information, some care is needed when 
interpreting the results of a /3 analysis. In the lower left-hand triangle of Table 5 we 
list the (3s derived by comparing the "sky" rmses of each pair of channels. To estimate 
combined-channel (3s we ignore the Ka4 channel, since we know that the "sky" rms for this 
channel is significantly higher than that estimated from the likelihood analyses, and average 
over all the other (3s. We find (3^3. ~ —0.2, (3q 2.2, and /^Ka+Q ~ 0.9. These should 
be compared to those derived in G95 from likelihood analyses using the fiat bandpower 
spectrum: (3^^ = 0.2^°'^, (3q = 1.7^J'g, and /^Ka+Q = 0.9^°'g. Our approximate estimates 
of the (3 central values reproduce the more accurate estimates of G95 to within cr/3, where 
the cr are the G95 error bars on the (3 values. This agreement gives some indication of 
the importance of off-diagonal noise and spatial correlations, as well as of the effect of 
oversampling. 

The upper right-hand triangle of Table 5 gives the (3s which should be found if the 
sky were described by a flat bandpower CMB spectrum. Averaging as before, but this 
time including the Ka4 channel, we flnd /^Ka ~ 0.5, (3q 0.4, and /^Ka+Q ~ 0.4. A similar 
analysis was carried through for all the other model CMB angular spectra used here, with 
the result that /^Ka ~ 0.5 — 0.8, (3q ~ 0.4 — 0.6, and /^Ka+Q ~ 0.4 — 0.7, with the steeper 
slopes corresponding to the spatially-flat A CDM models. That is, since the SP94 window 
functions couple spatial and frequency information one does not expect to flnd (3 = even 
for a pure CMB spectrum. Our approximate (3s agree with the values found by G95 using 
flat bandpower likelihood analyses to ~ cr/S for the Ka data and to ~ O.Scr for the Q and 
Ka -|- Q data. We therefore conclude that the G95 data could be even more consistent with 
a CMB spectrum than was indicated there. Since the SP94 Wi couple /- and i/-space, the 
only way to see how consistent the G95 (3 values are with a CMB anisotropy spectrum is to 
use a model for both the spatial (/) and frequency (u) dependence of the assumed non-CMB 
foregrounds. While not statistically signiflcant at even the Icr level it is noted that for all 
CMB spectra we consider (3q < /^Ka- 

As shown in the second last column of Table 4, the deduced Ka -|- Q data average 
absolute uncertainties are larger than those for the Ka data and only slightly smaller than 
those for the Q data. Provided that there is no signiflcant non-CMB contamination in 
one of the bands and that the absolute-pointing uncertainty can be ignored, these error 
bars only account for instrumental and atmospheric noise, which must integrate down, and 
sample variance due to the limited number of independent pixels of data, which will not 
integrate down without more, independent, sky coverage. 

To investigate the behaviour of the error bars, in Table 6 we list the bandtemperature 
central values and ET and HPD limits derived for each channel using a flat bandpower 
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spectrum. The Ql channel data has a detection for flat bandpower, as well as for some 
of the low-density open model CMB spectra. However, accounting for calibration and 
beamwidth uncertainties, the flat bandpower model —2(7 HPD 5Ti limit is ~ 3/iK, so it is 
not a robust detection. The Ka2 channel does not have a detection, independent of the 
assumed CMB anisotropy spectrum. 

We may approximately model the total average absolute error bars for the combined 
Ka and Q data sets of Table 4, and for the individual channel data sets of Table 6, using 

O"tot.,com. = ^CTSV^ + CT^'^ / Nj 
O"tot.,ind. = \/o^V^~+~0^, (23) 

where ctsv is the uncertainty due to sample variance, ctn that due to intrinsic noise, and 
Ni is the number of channels that contribute to the corresponding combined data set. We 
emphasize that the model of eqs. (23) is approximate, since it neglects correlations between 
channels. 

To determine cr^^ -^^^ , where A is either Ka or Q, we use cr^^ -^^^ = [Y^ilj^iiu^^Y / Nj]^^"^ , 
where the are the individual-channel HPD STi average absolute error bars of Table 6. 
With robust detections only, crj^t^ind. = 14.6 /iK (iVj = 3) and (j2,t.,ind. = 12.1 /iK [Ni = 2). 
With all channels, af^l^.^d. = 13-8 /iK [Ni = 4) and (j2,t.,ind. = 12-3 /iK [Ni = 3). 

Using these numbers and the average HPD 5Ti error bars of Table 4, 10.1 /iK for Ka 
and 11.7 /iK for Q, to explain the behaviour of the Ka and Q data error bars the sample 
variance and intrinsic noise contributions to the 5Ti error bars need to be 

(jf^ = 6.8(8.5) /iK, (jj^^ = 13(11) /iK, 

a2y = 11(11) /iK, = 4.4(4.6) /iK. (24) 

Here the term in parentheses on the right hand side of each equation is that determined 
using all channels, and the flrst terms are those determined using just the robust detection 
channels. This approximate estimate suggests that the Q data total error bars are 
dominated by sample variance, the Ka data total error bars are dominated by intrinsic 
noise, and that the Q data is intrinsically less noisy than the Ka data. The behaviour of 
the error bars (in Table 4) when the Ka and Q data are combined to form the Ka -|- Q data 
set then does not seem very surprising. 

It is possible to simply estimate the sample variance error bars, though a more accurate 
estimate would make use of simulations, and would be somewhat model dependent. We 
flrst illustrate the method using the SK93 overlap part of the SK94 experiment (N95), since 
N95 have performed the requisite simulations and we can use their result to judge how 
accurate the estimate is. 
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Most of the weight of the SK93 overlap data comes from the Ka-band observations (L. 
Page, private communication 1995), which were performed with a beamwidth Cp^jji^ = 1.4°. 
The beam centre traces out a circle of radius 4.4° centered at the NCP, so the experiment 
has A^pix = {2tt x 4.4°)/1.4° = 20 pixels (the data is binned into 24 bins, so the SK93 overlap 
experiment somewhat oversamples). The sample variance error is l/y'2iVpi^ = 16%. 
From the Ka93 entry in Table 1 of N95, Arms = 37 /iK, so the approximate estimate of the 
sample variance error bar is 5.9 /iK, which should be compared to the N95 result of 6.5 /iK 
(last paragraph on p. L71; note that this is the sample variance contribution to the ST^^s 
error bars). The simple analytical estimate of sample variance needs to be increased by 
10% to reproduce the N95 result. 

Although there are a number of differences between the SK93 overlap and SP94 
experiments, the analytical method to estimate sample variance should be fairly accurate 
for SP94. At the central Ka-band and Q-band SP94 frequencies, the beamwidths are 
4wHM = 1-47° and 4whm = 1-11°, so for a 20° scan A^g^^ = 13.6 and N^^J = 18.0, so 
SP94 is more oversampled than the SK93 overlap. Accounting for the 10% correction, the 
SP94 Ka sample variance is 21.1%, and the Q sample variance is 18.4%. Using the Ka and 
Q 5Ti HPD central values of Table 4, 30.4 /iK and 41.4 /iK respectively, and eqs. (23), we 
find for the sample variance and intrinsic noise contributions to the SP94 5Ti average error 
bars. 



(7. 
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6.4(6.4) /iK, (7^^^ = 13(12) /iK 



(j^v = 7.6(7.6) /iK, = 9.4(9.7) /iK, (25) 

which should be compared to the estimates of eqs. (24). Again, the first term on the right 
hand side of each equation is that determined using just the robust detection channels, and 
the terms in parentheses use data from all channels. This approximate estimate of the Ka 
sample variance and intrinsic noise is close to what is required (eqs. [24]) to explain the 
behaviour of the Ka data error bars. However, the approximate estimate of the Q sample 
variance indicates that it is not as important (relative to the Q noise) as what is required 
(eqs. [24]) to explain the behaviour of the Q data error bars. Given the small sky coverage 
of SP94, this discrepancy might not be statistically significant. More data, with greater 
sky coverage, is needed to resolve this issue. However, it is clear that sample variance is a 
significant contributor to the SP94 error bars. 

In Table 7 we give bandtemperature values for the combined Ka, Q and Ka -|- Q 
data subsets corrected for beamwidth and calibration uncertainties in two different ways. 
The corresponding uncorrected numbers are given in Table 4. As mentioned previously, 
we have ignored the absolute-pointing uncertainty when combining data from the Ka and 
Q scans, and for this Ka -|- Q data set analysis we have also assumed that the Ka and Q 
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beamwidth uncertainties shift together. These effects may approximately compensate, but 
we emphasize that the results from the Ka + Q data set analysis are less robust than those 
from either the Ka or Q data subset analyses. 

The approximate prescription generally used in the CMB anisotropy field to account 
for these uncertainties results in total la error bar ranges that are not very different from 
the total 1(7 HPD error bar ranges derived here (rows labelled JOG vs. rows labelled HPD 
of Table 7). Note, however, that this is because the upper HPD error bars derived here are 
somewhat smaller and the lower HPD error bars are somewhat larger than the error bars 
previously used to represent the SP94 combined Ka and Q results (RSBG; GRS). Given 
our understanding of the numerical uncertainties of the computations here, we believe that 
these differences are numerically significant. If a similar result holds for a number of other 
small-scale CMB spatial anisotropy experiment data analysis results, then this would raise 
the lower values of reduced-^^ found by GRS and probably also tighten the range in the 
values for a given model. 

Tables 8 and 9 show the effects of the beamwidth- and calibration-uncertainty 
corrections on the deduced Qrms-PS central values and limits for some of the SP94 data 
subsets. The corresponding likelihood functions for the selected models are shown in Figs. 5 
and 6. 

The beamwidth-uncertainty correction has a bigger effect for those data sets that probe 
a region of /-space where the model C; spectra changes rapidly with /. This is because 
in this case the contributions from the iblcr beamwidth analyses no longer approximately 
compensate for each other. This effect is clearly evident in the numbers of Tables 8 and 9. 

As discussed in the previous section, the calibration uncertainty correction broadens 
and skews the probability density distribution functions toward higher values of Qrms-PS, 
resulting in a small increase in the central Qrms-PS values and a larger increase in the la 
limit range. These effects are also clearly evident in the numbers of Tables 8 and 9. 

Tables 10 and 11 show the effects of the beamwidth- and calibration-uncertainty 
corrections on the maximum and marginalized values of the probability density distribution 
functions for some of the SP94 data subsets. Note that the Ka2 data set does not have 
a detection. From Table 11, we see that these corrections do not significantly affect the 
ordering of models with respect to the marginal probability distribution function value. 

It is interesting that at a given beamwidth- and calibration-uncertainty treatment, 
there is a greater range in marginal probability values for the Ka data then there is for the 
Ka -|- Q data (see Table 11); this larger range is also evident for the Q data (see Table 
16). In all cases, among all the models considered here, low-density Oq ~ 0.1 — 0.2 open 
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models (models 01 - 04) are favoured by the SP94 data. If we assume that the marginal 
probability distribution (the probability density distribution for just the cosmological 
parameters Oq, ^, and 0^) is gaussian, then relative to a model with a renormalized 
marginal value of unity, a model Icr away has a value = 0.61, and one 1.5(7 away has a value 
= 0.32. Clearly, the SP94 data do distinguish between the model CMB anisotropy shapes, 
but not at a very significant level. 

Tables 12 - 16 summarize our SP94 data set analyses. 

Table 14 lists the central values and Icr ranges of band temperature for the individual- 
channel data sets. For a given data set the bandtemperature values vary by ~ — 10% from 
model to model, with a typical range of ^ 5%. Given the other uncertainties, an analysis 
based on the fiat bandpower CMB spectrum is still quite useful for deriving approximate 
constraints on models, but it is not ideal. This is because a fiat bandpower analysis cannot 
provide information about the relative probability of the model CMB anisotropy shapes. 
Since the combined Ka, Q, and (full) Ka + Q data subset analyses are based on data from 
a combination of windows, it is not possible to accurately compute the corresponding 5Ti 
central values and limits, except for the fiat bandpower model (and these are given in Table 

7). 

Note that the variation in 5Ti with beamwidth is approximately — 4 percentage points 
less than the variation in Qrms-PS with beamwidth which, depending on model and data 
subset, can be ~ 15% between the —Icr and -|-lcr beamwidths. Hence the variation in 5Ti 
with beamwidth is still quite significant. This is particularly true for the fiat bandpower 
spectrum, where 5Ti varies as much as Qrms-PS, and so does not take out any of the effect 
of Wi varying with beamwidth. 

Focussing on just the SP94 Qrms-PS central values and limits from the combined data 
sets (last column of Table 12 and last two columns of Table 13), one notices that for each 
model the Q-data Qrms-PS values are somewhat higher than those for the Ka data. However, 
this is only at approximately Icr of either the Ka or Q data error bars, and is not significant. 
Comparing the corresponding 5Ti values for the fiat bandpower model (fourth, seventh, and 
tenth row of Table 7) one again notices this effect, but we emphasize that given the errors 
this is not significant. Furthermore, the Ka + Q data set values are consistent with those 
from the Ka and Q data subsets, although the Ka + Q data values are much closer to the 
Q data values. 

The individual-channel 5Ti central values and ranges (Tables 14) are approximately 
model independent for each data set, and so provide a convenient, approximate summary of 
the SP94 observational results. Focussing on a given model and recalling that the frequency 



- 24 - 



dependent SP94 window function couples different channels to different regions of /-space, 
one sees that, in a given line in Table f4, the SP94 STis seem to mildly drop or remain 
constant from / ~ 40 to / ~ 50. They then seem to mildly rise from / ~ 50 to / ~ 60, or at 
least stay approximately constant but at a higher value than for the 40 ^ / ^ 50 range. We 
emphasize that the precise numerical / values are model dependent, and that the quoted 
numbers correspond approximately to those for model Of. We also note that the higher 
values at larger / are from the Q data, so it is prudent to bear in mind the earlier discussions 
about (3 and sample variance. This effect may be seen more clearly in Figs. 7, and is related 
to the effect mentioned in the previous paragraph; again, this is probably not significant. It 
is interesting that the mildly falling or fiat 5Ti values at lower / are qualitatively consistent 
with what is seen in the SK94 and SK95 experiments (Netterfield 1995; Netterfield et al. 
1997), and that the "dip" or "break" at / ~ 50 — 60 is also qualitatively consistent with 
what is seen in the SK95 experiment (Netterfield et al. 1997).^ It would be of interest to 
determine whether or not the SP94 and SK95 "dips" are quantitatively consistent. We 
emphasize that here we have assumed that the SP94 data is purely CMB anisotropy. The 
"dip" however might be a consequence of as yet undiscovered small systematic effects, 
foreground contamination, statistical noise, neglect of the absolute pointing uncertainty, 
sample variance, or just a result of the CMB anisotropy not being well-described by the 
models we have used here. Even if there is a "dip" in the CMB anisotropy, it is not very 
significant in the SP94 data. More data will be needed to resolve these issues. 

The marginal probability distribution values for the combined SP94 data sets (last 
three columns of Table 16) are similar, although there is a larger range for the Ka and Q 
data sets than for the Ka -|- Q data set. CMB anisotropy shapes that are favoured by the 
Ka data are also favoured by the Q data and hence by the Ka -|- Q data. Those that are 
not favoured by the Q data are also not favoured by the Ka data and hence neither by the 
Ka -|- Q data. 

The Ka -|- Q data, accounting for beamwidth- and calibration-uncertainty, favour (last 
column of Table 16) an open Oq = O.f model (model Of), among the models we consider 
here. We do not have the resolution in model-parameter space to construct the complete 
marginal probability distribution function (the probability density distribution function for 
just the cosmological parameters), and so cannot properly judge the significance of the 
values we have computed at isolated points in model-parameter space. To get some idea of 
the significance of these values we may assume that the marginal distribution is not far from 
a gaussian. As mentioned above, for a gaussian marginal distribution, a model fcr away 



^Note that in Figs. 7 - 9 of Netterfield et al. (1997) the SK data are plotted at /g of the Wi, and not at 
of the ((^Trms^); as done in Fig. 7 here. 
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from the most favoured low-density open model CMB anisotropy shape has a marginal 
value of 0.61, and a model with marginal value 0.32 is 1.5(7 away from the most favoured 
low-density open model. If the marginal probability distribution function is narrower than 
a gaussian the above cr-values are overconservative, and if it is wider than a gaussian they 
are overoptimistic for the purpose of ruling out models. Under the gaussian assumption, 
the difference between the most favoured low-density open model and the least favoured 
model, which is always a flat-A one, is ~ l.Gcr for Ka -|- Q and ~ l.Scr for Ka or Q, which 
are not very significant. 

We conclude that the SP94 data are most consistent with the CMB anisotropy shape 
in low-density open CDM models with Oq ~ 0.1 - 0.3 and 0.4 with larger h and smaller 
Ob, and with the fiat bandpower shape, at least among the models we consider here. Since 
the fiducial CDM and fiat-A models have CMB anisotropy shapes that are always more 
than Icr away, under the gaussian assumption, from the most favoured low-density open 
model, the SP94 data do not favour these models. This is especially true for old (to ^ 15 - 
16 Gyr), large baryon density (0^ ^ 0.0175/i~^), low-density (Oq ~ 0.2 - 0.4), fiat-A ones. 
These results are mostly consistent with the qualitative ones of RSBG, and the quantitative 
ones of GRS. We emphasize, however, that under the gaussian marginal assumption the 
SP94 data do not rule out any of the models we consider here at the 2cr level, on the basis 
of their CMB anisotropy shape alone. 

Finally, we note that the model normalizations deduced from the two-year DMR data 
and from the SP94 data sets are mostly consistent although, as mentioned above, the Ka 
data prefer a normalization ~ Icr below that favoured by the Q data. The SP94 Q and 
Ka -|- Q data favour a slightly higher normalization for the Oq = 0.1 open model than 
do the two-year DMR data, and the SP94 Ka and Ka -|- Q data favour a slightly lower 
normalization of the older (smaller /i), higher 0^, low-density fiat-A models than do the 
two-year DMR data. These results are consistent with those of GRS. It will be of interest 
to see if the Oq = 0.1 open model can be significantly ruled out on the basis of CMB 
anisotropy data alone when the new, slightly lower, four-year DMR normalization for this 
model is combined with our results here. 

4. Conclusion 

We have developed general methods to account for various sources of experimental 
and observational uncertainty in the likelihood analysis of observational data from a CMB 
spatial anisotropy experiment. 
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We have accounted for beamwidth- and calibration-uncertainty in likelihood analyses 
of the SP94 observational data that make use of theoretical CMB spatial anisotropy spectra 
in a variety of observationally-motivated, open and spatially-flat A, CDM cosmogonies. 

Absolute-radiometer-calibration and beamwidth uncertainties are the largest known 
sources of uncertainty for the Ka- and Q-scan data from the SP94 experiment that were 
not previously explicitly accounted for. In our analysis we have not explicitly accounted 
for a number of other known but very much smaller uncertainties, including: the ellipticity 
and possible mild nongaussianity of the beams; the relative beamwidth uncertainty of the 
individual-channel beams; the chop-amplitude uncertainty; the effect of the main beam 
efficiency; the relative-radiometer-calibration uncertainty of the individual channels; the 
uncertainty in the measured passband frequencies; and the relative-pointing uncertainty. 
Collectively, these are likely to contribute to the error bars at the percentage-point level. In 
addition, when combining the Ka- and Q-scan data we have ignored the absolute-pointing 
uncertainty. We have also analyzed the Ka -|- Q data set assuming that the Ka and Q 
beamwidth uncertainties shift together. These two effects, quite likely, compensate, and our 
Ka -|- Q data results are not grossly incorrect; however, the Ka and Q data subset results 
are more robust. 

In our analyses of the SP94 data we have assumed that it is purely CMB spatial 
anisotropy; justification for this assumption may be found in G95 and Gundersen (1995). 

The combined Ka, Q, and (full) Ka -|- Q data subsets show evidence for at least 
3(7 detections of anisotropy. Six of the seven individual-channel "sky" rms values agree 
with the rms values derived from individual-channel fiat bandpower analyses to within 
25%. Predicted fiat bandpower (3 values agree with those found from fiat bandpower 
likelihood analyses by G95 to within a third (Ka) and 0.8 (Q and Ka -|- Q) of the G95 
error bars. Hence the G95 data could be even more consistent with a CMB i/-spectrum 
than was indicated there. However, since the SP94 Wi{i') couple /- and i/-space, the only 
way of establishing from a /3-analysis how inconsistent the SP94 data is with non-CMB 
contamination would be to use a model for both the /- and i/-dependence of the suspected 
non-CMB foreground in likelihood analyses of the data. Only if a multiple-frequency CMB 
anisotropy experiment has individual-channel window shapes (in /) that are identical does 
(3^0 indicate that the data are inconsistent with a CMB i/-spectrum. The behaviour of 
the Ka data error bars is consistent with what is expected from our approximate estimate 
of sample variance; the behaviour of the Q data error bars is not as consistent. For all the 
model CMB spatial anisotropy spectra we consider here, the Q data prefer a somewhat 
higher normalization than do the Ka data, but this is at only ~ Icr of either the Ka- or 
Q-data error bars. To see if the differences between the Ka and Q data are significant will 
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require more data. Until this is accomplished it would be prudent to not draw conclusions 
using only the full Ka + Q analysis, but to also consider the implications of both the Ka 
and Q analyses. 

The HPD bandtemperature estimates of Table 7 provide an approximate summary of 
the SP94 data, and can be used to approximately constrain other models of CMB spatial 
anisotropy, although likelihood analyses, like those done here, are needed if one wishes to 
utilize all the information in the SP94 data. Furthermore, as discussed above, the deduced 
bandtemperature values typically vary by ~ 5% from model to model. 

While the marginal probability distribution function values show that the SP94 data 
are clearly sensitive to spatial structure on the sky, they do not significantly distinguish 
between models, although they do favour low-density open models over older, high 0^, 
low-density, fiat- A models. No model considered here is ruled out at the 2cr level by the 
SP94 data alone, at least in the gaussian marginal probability distribution approximation. 
The SP94 model normalizations are mostly consistent with those deduced from the 
two-year DMR data. These results are mostly compatible with what one concludes from a 
goodness-of-fit analysis of all presently available CMB anisotropy detection data (GRS). 

The major contributors to the SP94 Ka- and Q-scan data analyses error bars are 
(alphabetically): (1) beamwidth uncertainty; (2) calibration uncertainty; (3) intrinsic noise 
uncertainty; and, (4) sample variance uncertainty. More sky coverage and less oversampling 
will reduce the sample variance contribution, and longer integration will reduce the intrinsic 
noise contribution. Measured beamwidth and calibration distribution functions should 
help reduce the beamwidth- and calibration-uncertainty contributions, unless it turns out 
that our gaussian approximation is actually narrower than the real distribution function. 
It should be possible, in the next few years, to acquire data with error bars significantly 
smaller than the present ones (~ 30 — 36%). 

If measured distribution functions for the parameters of CMB experiments, at least for 
those that are known to have relatively large uncertainty, become available, the methods 
we have developed here will allow for a more robust estimate of the sky signal detected in 
CMB anisotropy experiments, even for presently available CMB anisotropy data. 

We acknowledge helpful discussions with D. Bond, K. Gorski, J. Peebles and the 
referee, E. Wright, and are especially indebted to L. Page. KG acknowledges support 
from NASA grants NAGW-4623 and 4870. BR acknowledges support from NSF grant 
EPS-9550487 and matching support from the state of Kansas. 
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Table 1: Numerical Values for the Individual-Channel Zero-Lag Window Function 
Parameters^ 





Chan.: 


Ivai 


Ivaz 


Ivao 


TZ^ /I 


(41 


(4Z 


(46 


Beamwidth 


i/(GHz): 


27.25 


29.75 


32.25 


34.75 


39.15 


41.45 


43.75 




/g-0.5 


34 


35 


37 


38 


40 


41 


42 






55.2 


58.0 


60.5 


62.8 


67.1 


68.8 


70.4 


-1(7 




62 


65 


67 


70 


73 


75 


76 




/g-0.5 


95 


100 


104 


107 


114 


116 


118 






1.04 


1.09 


1.14 


1.17 


1.24 


1.27 


1.29 




/g-0.5 


33 


34 


36 


37 


39 


40 


40 




/e 


53.4 


56.1 


58.7 


61.0 


64.4 


66.1 


67.8 


Nom. 


/m 


60 


63 


65 


68 


71 


73 


74 




/g-0.5 


92 


97 


101 


105 


110 


112 


115 






1.01 


1.06 


1.10 


1.14 


1.20 


1.23 


1.26 




/g-0.5 


32 


33 


35 


36 


38 


38 


39 




/e 


51.6 


54.4 


56.9 


59.2 


61.9 


63.7 


65.3 


+ la 


/m 


58 


61 


64 


66 


69 


70 


72 




/g-0.5 


89 


94 


98 


102 


106 


109 


111 






0.997 


1.03 


1.07 


1.11 


1.16 


1.19 


1.22 



^Computed according to the prescriptions above and in eq. (4). 
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Table 2: Numerical Values for Parameters Characterizing the Shape of [ST^ 



Window: 






+ lcr Kal 






-la Q3 




# 




/g-0.5 




/g-0.5 


/g-0.5 




/g-0.5 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


01 


(0.1, 0.75, 0.0125) 


7 


38 


74 


14 


55 


102 


02 


(0.2, 0.65, 0.0175) 


14 


43 


79 


22 


63 


108 


03 


(0.2, 0.70, 0.0125) 


13 


42 


78 


21 


61 


107 


04 


(0.2, 0.75, 0.0075) 


13 


41 


77 


20 


59 


104 


05 


(0.3, 0.60, 0.0175) 


17 


47 


81 


27 


66 


111 


06 


(0.3, 0.65, 0.0125) 


16 


46 


80 


25 


64 


108 


07 


(0.3, 0.70, 0.0075) 


16 


44 


78 


24 


61 


106 


08 


(0.4, 0.60, 0.0175) 


19 


48 


82 


29 


67 


111 


09 


(0.4, 0.65, 0.0125) 


18 


47 


80 


27 


64 


109 


OlO 


(0.4, 0.70, 0.0075) 


18 


45 


79 


25 


62 


106 


Oil 


(0.5, 0.55, 0.0175) 


21 


50 


83 


31 


68 


113 


012 


(0.5, 0.60, 0.0125) 


20 


48 


81 


29 


66 


110 


013 


(0.5, 0.65, 0.0075) 


19 


46 


79 


27 


63 


107 


014 


(1.0, 0.50, 0.0125) 


22 


49 


83 


31 


68 


115 


Al 


(0.1, 0.90, 0.0125) 


28 


57 


90 


40 


76 


121 


A2 


(0.2, 0.80, 0.0075) 


24 


53 


87 


35 


73 


120 


A3 


(0.2, 0.75, 0.0125) 


26 


55 


88 


37 


75 


121 


A4 


(0.2, 0.70, 0.0175) 


27 


56 


89 


39 


76 


122 


A5 


(0.3, 0.70, 0.0075) 


23 


52 


86 


34 


72 


119 


A6 


(0.3, 0.65, 0.0125) 


25 


54 


88 


36 


74 


120 


A7 


(0.3, 0.60, 0.0175) 


27 


56 


89 


38 


75 


121 


A8 


(0.4, 0.65, 0.0075) 


23 


51 


85 


32 


70 


118 


A9 


(0.4, 0.60, 0.0125) 


24 


53 


87 


35 


73 


119 


AlO 


(0.4, 0.55, 0.0175) 


26 


55 


88 


37 


75 


121 


All 


(0.5, 0.60, 0.0125) 


23 


51 


86 


33 


71 


118 


A12 


(1.0, 0.50, 0.0125) 


22 


49 


83 


31 


68 


115 


Flat 




16 


41 


73 


21 


54 


96 


Wi 




32^ 


51.6*^ 


89^ 


42a 


70.4*^ 


118^ 



^/g-0.5 for the window, computed according to the prescription above eq. (4). 
^le for the window. 
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Table 3: Numerical Values for Rms Temperature Anisotropies' 



Channel: 


Kal 


Ka2 


Ka3 


Ka4 


Ql 


Q2 


Q3 


"Sky"b 


39 


18 


36 


51 


33 


47 


42 


FBP, Nom. Wi"" 


51 


17 


37 


33 


34 


49 


52 



'^Estimated from the data of Fig. 4, as discussed in §2. 

■^Converted to rms (using eq. [7]) from the results of the likelihood analysis for the flat bandpower (FBP) 
angular spectrum (eq. [3]), and ignoring beamwidth and calibration uncertainties. 



Table 4: Numerical Values for Bandtemperature^from Likelihood Analyses Assuming a Flat 
Bandpower Spectrum and Ignoring Beamwidth and Calibration Uncertainties 



Data Set 


Analysis 


-1(7 


Peak 


+ 1(7 


Ave. Abs. Err.'' 


Ave. Fra. Frr.^ 




095*^ 


24 


29 


43 


±10 


±33% 


Ka 


Here" 


25 


30 


46 


±11 


±35% 




Here** 


22 


30 


42 


±10 


±33% 




095*^ 


33 


40 


55 


±11 


±28% 


Q 


Here*^ 


35 


41 


59 


±12 


±29% 




Here** 


32 


41 


55 


±12 


±28% 




095*^ 


30 


36 


50 


±10 


±27% 


Ka + Q 


Here*^ 


33 


39 


56 


±11 


±29% 




Here** 


30 


39 


52 


±11 


±27% 



""STi (eq. [7]) in //K. 
•^ET (16.0% & 84.0%) prescription. 
=ET (15.87% & 84.13%) prescription. 
■^HPD (68.27%) prescription. 
Average absolute error in //K. 

^Average fractional error, as a fraction of the central value. 
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Table 5: Numerical Values for /3, from the "Sky" Rms^ and from a Theoretical Flat 
Bandpower CMB Angular Spectrum*^ 



Channel 


Kal 


Ka2 


Ka3 


Ka4 


Ql 


Q2 


Q3 


Kal 




0.55 


0.53 


0.52 


0.48 


0.47 


0.46 


Ka2 


-8.7 




0.51 


0.50 


0.46 


0.45 


0.44 


Ka3 


-0.47 


8.6 




0.48 


0.44 


0.43 


0.42 


Ka4 


1.1 


6.7 


4.6 




0.42 


0.41 


0.40 


Ql 


-0.49 


2.1 


-0.51 


-3.7 




0.40 


0.39 


Q2 


0.43 


2.9 


1.0 


-0.48 


6.3 




0.38 


Q3 


0.15 


2.2 


0.49 


-0.85 


2.2 


-2.1 





^Lower left-hand triangle. Estimated from the ratio of the signal in the channel of the first row, to the signal 

in the channel of the first column (eq. [21]). 
'^Upper right-hand triangle. Estimated from the ratio of the signal in the channel of the first column, to the 

signal in the channel of the first row (eq. [22]). 
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Table 6: Numerical Values for Bandtemperature^from Likelihood Analyses Assuming a Flat 
Bandpower Spectrum, for the Individual-Channel Data Sets, and Ignoring Beamwidth and 
Calibration Uncertainties 



Data Set 


Prescript. 


-1(7 


Peak 


+ la 


Ave. Abs. 


Err.*^ Ave. Era. Err.*^ 


Kal 


ET 


40 


51 


79 


±20 


±39% 




HPD 


35 


51 


73 


±19 


±37% 


Ka2 


ET 


7.7 


16 


30 


±11 


±69% 




HPD 


5.4 


16 


28 


±11 


±67% 


Ka3 


ET 


27 


34 


53 


±13 


±40% 




HPD 


23 


34 


49 


±13 


±38% 


Ka4 


ET 


23 


29 


45 


±11 


±39% 




HPD 


20 


29 


42 


±11 


±37% 


Ql 


ET 


19 


28 


45 


±13 


±47% 




HPD 


17 


28 


42 


±13 


±46% 


Q2 


ET 


33 


39 


58 


±12 


±32% 




HPD 


30 


39 


53 


±12 


±30% 


Q3 


ET 


35 


42 


60 


±13 


±31% 




HPD 


31 


42 


56 


±12 


±29% 



""STi (eq. [7]) in //K. ET limits are 15.87% & 84.13% and HPD limits are 68.27%. 
'^Average absolute error in //K. 

■^Average fractional error, as a fraction of the central value. 
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Table 7: Numerical Values for Bandtemperature^from Likelihood Analyses Assuming a Flat 
Bandpower Spectrum and Accounting for Beamwidtli and Calibration Uncertainties 



Data Set 


Prescript. 


-1(7 


Peak 


+ la 


Ave. Abs. 


Err.** Ave. Era. Err.'' 




JOG*^ 


22 


29 


44 


±11 


±37% 


Ka 


ET" 


25 


31 


48 


±12 


±38% 




HPD-^ 


22 


31 


43 


±11 


±36% 




JOG*^ 


31 


40 


57 


±13 


±32% 


Q 




35 


42 


61 


±14 


±32% 




HPD-^ 


31 


42 


57 


±13 


±31% 




JOG*^ 


28 


36 


51 


±11 


±31% 


Ka + Q 




33 


39 


58 


±13 


±32% 




HPD-^ 


30 


39 


53 


±12 


±30% 



""STi (eq. [7]) in //K. 

'^JOG prescription to account for calibration and beamwidth uncertainties by adding 15% of the central value 
in quadrature to the G95 ET (16.0% & 84.0%) prescription error bars. 

■^ET (15.87% & 84.13%) and HPD (68.27%) limits from the beamwidth- and calibration-uncertainty corrected 

probability density distribution functions. 
■^Average absolute error in /iK. 

Average fractional error, as a fraction of the central value. 
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Table 8: Numerical Values for Qrms-PS (in /^K) from the Ka2 and Q2 Data' 









Data 


Set: 




Ka2 






Q2 




Model 




(Oo, ^Bh^) 


Tt b 

Bare 


Beam. 


B.+C. 


Tt b 

Bare 


T* r 

r>eam. 


B.+C. 


01 


(0 


1 





75 





0125) 


17 . 


1 T 43 

17 .: 


17 . 


Q 48 
00 27 


36 4« 


0/^ 49 


02 


(0 


2 





65 





0175) 


16 .! 


^ r A^ 

16 .! 


^ r A'^ 

16 .: 


Q c 47 
•'>'^ 26 


Q C 47 
'J'^ 26 


Q r 48 


03 


(0 


2 





70 





0125) 


1 T 42 

17 . 


1 T 42 

17 . 


-1 1-! AA 

17 . 


0/^ 48 

36 4« 


0/^ 48 

36 4« 


0/^ 49 


04 


(0 


2 





75 





0075) 


1 T 43 

17 .: 


1 T 43 

17 .: 


-1 ^7 45 

17 .? 


07 49 

0/ 27 


07 49 
0/ 27 


0751 

6t 27 


05 


(0 


3 





60 





0175) 


14;.: 


14:.: 


1 /1 37 
14:.. 


on 40 
30 23 


or\ 41 
30 23 


on 42 

30^^ 


06 


(0 


3 





65 





0125) 


14:. 


1 /I 36 

14:. 


1 /I 38 

14:. 


01 41 
"Jl^ 23 


01 42 
•'>^ 23 


01 43 

•'>^ 23 


07 


(0 


3 





70 





0075) 


1 r 37 

15:.! 


1 r 38 

15:. 


1 r 39 

15:. 


00 43 
OZ 24 


QO 43 
OZ 24 


Q r) 44 

0Z24 


08 


(0 


4 





60 





0175) 


1 ri 29 

12 . 


1 ri 29 

12 . 


1 ri 30 

12:. 


ri r 33 

253^ 


f; 34 
ZO 19 


ri r 35 

25^5 


09 


(0 


4 





65 





0125) 


1 ri 30 


1 ri 30 


1 ri 31 

12:.: 


2635 


2635 


26 ?6 


OlO 


(0 


4 





70 





0075) 


1 ri 31 

12:.: 


1 ri 31 

12:.: 


1233 


27^15 


27^15 


07 37 
^ ' 20 


Oil 


(0 


5 





55 





0175) 


r> ri 23 

9.2 _ 


r> ri 23 

9.2 _ 


n 24 

9.3 _ 


rin 27 

2027 


rin 27 

2027 


on 28 

20 


012 


(0 


5 





60 





0125) 


9.6 _ 


f\ n 2A 

9.6 _ 


r> T 25 

9.7 _ 


rii 28 

21 2« 


rii 28 

21 2« 


rii 29 

21 16 


013 


(0 


5 





65 





0075) 


1 n 25 

10 .: 


1 n 25 

10 .: 


1 n 27 

10 .! 


22 2« 


00 29 

22 2« 


22^0 


014 


(1 








50 





0125) 


o o 21 

8.3 _ 


21 
8.3 _ 


/I 22 

8.4 _ 


1 24 
1824 


1 24 
1824 


1 25 
1825 


Al 


(0 


1 





90 





0125) 


n o 24 

9.8 _ 


n 24 

9.8 _ 


r> r> 25 

9.9 .: 


22 2« 


00 29 

22 2« 


22^0 


A2 


(0 


2 





80 





0075) 


9.6 _ 


f\ n 2A 

9.6 _ 


r> T 25 

9.7 _ 


rii 28 

21 2« 


rii 28 

21 2« 


rii 29 

21 16 


A ■ ■> 

A3 


(0 


2 





75 





0125) 


9.0 _ 


r\ n 99 


r\ 1 93 

9.1 . 


r*r\ 2fi 
2026 


r*r\ 27 

2027 


r*r\ 27 
2027 


A A 




2 





70 





0175) 


o.D ... 


O.D ... 


722 
0. t ... 


1 Q25 
19 14 


1 Q25 
19 14 


1 Q 26 
19 14 


A5 


(0 


3 





70 





0075) 


8.922 


8.922 


9.123 


20?^ 


2026 
15 


2027 
15 


A6 


(0 


3 





65 





0125) 


8.421 


8.421 


8.522 


18?^ 


19?^ 


19?^ 


A7 


(0 


3 





60 





0175) 


8.020 


8.020 


8.021 


18 f3 


18^3 


18?^ 


A8 


(0 


4 





65 





0075) 


8.621 


8.622 


8.722 


19?^ 


19?^ 


19?^ 


A9 


(0 


4 





60 





0125) 


8.120 


8.120 


8.221 


18^3 


18^3 


18?^ 


AlO 


(0 


4 





55 





0175) 


7.71.9 


7.71.9 


7.720 




17f3 


17?^ 


All 


(0 


5 





60 





0125) 


8.120 


8.120 


8.221 


18^3 


18^3 


18?^ 


A12 


(1 








50 





0125) 


8.321 


8.321 


8.422 


18^4 


18^4 


18?^ 


Flat 














1230 


1230 


1232 


Of; 34 
ZO 19 


2635 


2635 



^ The primary entry in each of the last six columns is the value at which the probability density distribution 
function peaks. Ellipses as the lower entry in a vertical pair denotes a non-detection; the corresponding 
upper entry is the 2(T (97.72% ET) upper limit. For detections, the vertical pair of numbers are the ±1(T 
(68.27% HPD) upper and lower limits. 

'^Computed using eq. (8). Ignores beamwidth and calibration uncertainties. 

■^Accounts for beamwidth uncertainty using eq. (12). Ignores calibration uncertainty. 

■^Computed using eqs. (12) and (16). Accounts for beamwidth and calibration uncertainties. 
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Table 9: Numerical Values for Qrms-PS (in /^K) from the Ka and Ka + Q Data' 









Data 


Set: 




Ka 






Ka + Q 




JViodel 




(Oo, h, neh^) 


TO b 

Bare 


T* r 

r>eam. 


B.+C. 


TO b 

Bare 


Beam. 


B.+C. 


01 
02 
03 
04 
05 


(0 
(0 
(0 
(0 
(0 


1 
2 
2 
2 
3 









75 
65 
70 
75 
60 









0125) 
0175) 
0125) 
0075) 
0175) 


^7 20 
97 37 
^ ' 20 
OO 38 
^^20 

28^^ 
23 


^7 20 
97 37 
^ ' 20 

00 38 
^^20 

28^^ 
23 


00 39 

^^20 

27 ^« 

00 39 

^^20 

00 40 

^^20 

23?? 


Q 46 
00 27 

34^^ 

K 45 
00 27 

0/^ 47 

36 S 

00 38 

23 


Q 46 
00 27 

. -) J 44 

34^^ 

Q r 46 

00 27 
0/^ 47 

36 S 

00 38 

23 


Q C 47 
'J'^ 26 
/I 46 

34^^ 

Q C 47 
'^'^ 26 
0/^ 48 

36^8 

00 39 
22 


06 
07 
08 
09 


(0 
(0 
(0 
(0 


3 
3 

4 
4 








65 
70 
60 
65 








0125) 
0075) 
0175) 
0125) 


ri ,1 33 

24 

c 34 
ZD 

1 r> 26 
rin 27 

20 2^ 


ri ,1 33 

24^7 

9 c 34 
ZD 

1 r> 26 

19 14 

rin 27 


24 ?4 

ri r 35 

25 ?| 

1 r> 27 

19?; 

on 28 

20?° 


on 39 

30 i 

Q 1 40 
Oi 24 

ri ,1 31 

24 ?i 

ri r 32 

25 ?2 


on 39 

30 i 

01 41 
Oi 24 

ri ,1 31 

24 ?i 

ri r 33 

25 ?3 


on 40 

30^^ 

01 42 
Oi 24 

/I 32 

24 ?2 

9 c 34 
ZD 19 


OlO 


(0 


4 





70 





0075) 


rii 28 

21 Is 


rii 29 
21^5 


rii 29 

21 ?5 


ri/^ 34 
2620 


ri/^ 34 

2620 


2620 


Oil 


(0 


5 





55 





0175) 


15 


15 fi 


16?? 


10 25 
15 


10 25 
15 


19 ?f 


012 


(0 


5 





60 





0125) 


12 


12 


16 23 
12 


20 26 


20 26 
ig 


20 ?I 


013 


(0 


5 





65 





0075) 


1723 
^ ' 12 


1723 
^ ' 12 


1724 
^ ' 12 


21 27 

^1- 16 


21 28 
^1- 16 


21 28 
^1- 16 


014 
Al 


(1 

(0 




1 






50 
90 






0125) 
0125) 


14 19 
10 
1 T 23 
17t2 


14 19 
10 
1 T 23 
17t2 


1420 
10 
1 T 23 

17?2 


1722 
1 ' 14 

rii 27 

21 ?6 


1723 
1 ' 14 
rii 27 

21 ?6 


1723 
1^ ' 13 
01 28 

21 ?6 


A2 


(0 


2 





80 





0075) 


16 


16 


16?^ 


20 26 


20 26 
ig 


20 ?I 


A3 


(0 


2 





75 





0125) 


15 fi 


15 fi 


15?? 


10 24 
15 


10 25 
15 


10 25 
15 


A A 


yj 


2 





70 





0175) 


1 A 20 

14ii 


1 A 20 

14ii 


1 f;2o 


1 S 23 
18 14 


1 S 23 
18 14 


1 8 24 
18 14 


A5 


(0 


3 





70 





0075) 


15- 


15 fi 


15?? 


19?^ 


19?^ 


19?^ 


A6 
A7 
A8 
A9 
AlO 


(0 
(0 
(0 
(0 
(0 


3 
3 

4 
4 
4 









65 
60 
65 
60 
55 









0125) 
0175) 
0075) 
0125) 
0175) 


10 
10 18 

i-J 10 

1420 
11 

14 19 

10 
10 18 
^■^ 10 


14 19 
10 

10 18 

i-J 10 

1420 
11 

14 19 

10 
10 18 
^■^ 10 


1420 
10 

14 19 
10 

1 520 

i-J 10 

141^0 
131' 


18 ?i 

17?^ 

18 ?i 
1722 

1^ ' 13 

16^1 

13 


18 ?i 

17?^ 

18 ?i 
1722 

1^ ' 13 

16^1 

13 


18 ?t 
17?^ 

18 ?t 
1723 

1^ ' 13 

16^1 

12 


All 
A12 
Flat 


(0 

(1 


5 







60 
50 






0125) 
0125) 


14 19 
10 

141^0 
20 ?I 


14 19 
10 

141^0 
20 ?I 


14 19 
10 

14?°o 
20 ?« 


1722 
^ ' 13 

1722 
i ' 14 

OK 33 
^"J 20 


1722 
^ ' 13 

1723 
i ' 14 

Of; 34 
^0 iQ 


1723 

^ ' 13 

1723 

^ ' 13 

Of; 34 
^0 19 



^Conventions are tlie same as for Table 8. 

'^Computed using eq. (8). Ignores beamwidth and calibration uncertainties. 
■^Accounts for beamwidth uncertainty using eq. (12). Ignores calibration uncertainty. 
■^Computed using eqs. (12) and (16). Accounts for beamwidth and calibration uncertainties. 
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Table 10: Renormalized Maximum Values of the Probability Density Distribution Functions' 



Data- 






Ka2 








02 






Ka 








Ka + Q 


Model 


Bare*^ 




B.+C.<* 






B.+C."* 






B.+C."* 


Bare*^ 






01 


1.0 


1.0 


1.0 


1 

X . \j 


1 

X . \j 


1 

X . \j 


1 

X . \j 


1 

X . \j 


1 

X . \j 


1.0 


1.0 


1 n 
i.U 


02 


n 
U 


00 

yy 


1.0 


1.0 


75 


75 


75 


81 


81 


81 


n 
U 


so 
oy 


n 
U 


so 
oy 


n so 
U.oy 


03 


n 
U 


00 

yy 


n 
U 


00 

yy 


n 
U 


00 

yy 


79 


79 


79 


84 


84 


85 


n 
U 


01 

yi 


n 
U 


01 

yi 


n 01 
u.y i 


04 


n 
U 


00 

yy 


n 
U 


00 

yy 


n 
U 


00 

yy 


0.84 


0.84 


0.85 


0.88 


0.88 


0.88 


n 
U 


00 

yz 


n 
U 


00 

yz 


n 00 
u.yz 


05 


n 
U 


OS 

yo 


n 
u 


OS 

yo 


n 
u 


OS 

yo 


0.63 


0.62 


0.62 


0.71 


0.71 


0.71 


n 
U 


so 
oZ 


n 
U 


so 
oZ 


n SI 

U.Oi 


06 


n 
U 


OS 

yo 


n 
U 


OS 

yo 


n 
U 


OS 

yo 


67 


68 


67 


75 


75 


75 


n 
U 


S/l 

o4 


n 
U 


S/l 

o4 


n S/l 

U.o4 


07 


n 
U 


OS 

yo 


n 
U 


OS 

yo 


n 
U 


OS 

yo 


74 


73 


74 


80 


80 


80 


n 
U 


Sfi 
oD 


n 
U 


Sfi 
oD 


n sfi 

U.oD 


08 


n 
U 


07 

y ( 


n 
U 


07 

y ( 


n 
U 


07 

y ( 


57 


57 


57 


66 


66 


66 


n 
U 


7S 
( o 


n 
U 


7S 
( o 


n 77 


09 


n 
U 


07 

y ( 


n 
U 


07 

y ( 


n 
U 


07 

y ( 


0.62 


0.62 


0.62 


0.71 


0.71 


0.71 


n 
U 


sn 
oU 


n 
U 


sn 
oU 


n sn 
u.ou 


OlO 


U 


Ofi 

yo 


n 
U 


07 

y ( 


n 
U 


07 

y ( 


69 


69 


69 


76 


76 


76 


U 


SQ 
00 


U 


SQ 
00 


n SQ 
U.oo 


01 1 


n 
U 


07 

y ( 


n 
U 


07 

y ( 


n 
U 


07 

y ( 


51 


51 


51 


61 


61 


61 


n 
U 


7/1 


n 
U 


7/1 
( 4 


n 7Q 

U. ( 


012 


n 
U 


Ofi 

yo 


n 
U 


07 

y ( 


n 
U 


07 

y ( 


57 


57 


57 


66 


66 


66 


n 
U 


77 


n 
U 


77 


n 77 


013 


n 
U 


Ofi 

yo 


n 


96 


n 


96 


64 


64 


64 


71 


71 


71 


n 
U 


sn 
oU 


n 
U 


sn 
oU 


n sn 

U.oU 


014 


n 


97 





97 





97 


0.49 


0.49 


0.49 


0.59 


0.58 


0.58 


n 


73 


n 


73 


n 70 
U. ( z 


Al 





96 





96 





96 


31 


31 


31 


41 


41 


41 





57 





57 


U.OD 


A2 





97 





97 





97 


39 


39 


39 


50 


49 

VJ . Tit./ 


49 

VJ . Tit./ 





65 





65 


65 


A3 





97 





97 





97 


0.35 


0.35 


0.35 


0.45 


0.45 


0.45 





61 





61 


0.60 


A4 





96 





96 





96 


0.32 


0.32 


0.32 


0.42 


0.42 


0.42 





58 





58 


0.57 


A5 





97 





97 





97 


0.42 


0.42 


0.42 


0.52 


0.52 


0.52 





68 





67 


0.67 


A6 





97 





97 





97 


0.37 


0.37 


0.37 


0.47 


0.47 


0.47 





63 





63 


0.62 


A7 





96 





96 





96 


0.33 


0.33 


0.33 


0.43 


0.43 


0.43 





59 





59 


0.59 


A8 





97 





97 





97 


0.45 


0.44 


0.44 


0.54 


0.54 


0.54 





70 





70 


0.69 


A9 





97 





97 





97 


0.39 


0.39 


0.39 


0.49 


0.49 


0.49 





65 





65 


0.64 


AlO 





96 





96 





96 


0.35 


0.35 


0.35 


0.45 


0.45 


0.45 





61 





61 


0.60 


All 





97 





97 





97 


0.43 


0.42 


0.42 


0.53 


0.53 


0.52 





68 





68 


0.67 


A12 





97 





97 





97 


0.49 


0.49 


0.49 


0.59 


0.58 


0.58 





73 





73 


0.72 


Flat 





94 





95 





95 


0.93 


0.94 


0.94 


0.92 


0.92 


0.93 





84 





84 


0.84 



^Renormalized such that it is unity for the model with the highest maximum value of the probability density 
distribution function, for the data set and the appropriate beamwidth- and calibration-uncertainty treatment. 
With the normalization set such that ^(Qrms-PS = 0) = 1, these highest maximum values of the likelihoods 
are roughly 2, 5 • 10^ 2 • 10^ and 5 • lO^^ for Ka2, Q2, Ka and Ka+Q, respectively. 

'^Computed using eq. (8). Ignores beamwidth and calibration uncertainties. 

■^Accounts for beamwidth uncertainty using eq. (12). Ignores calibration uncertainty. 

■^Computed using eqs. (12) and (16). Accounts for beamwidth and calibration uncertainties. 
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Table 11: Renormalized Marginal Values of the Probability Density Distribution Functions' 



Data- 






Ka2 












Q2 












Ka 










Ka + Q 


Model 


Bare*^ 




B.+C.<* 


Bare*^ 




B.+C.<* 


Bare*^ 




B.+C.<* 


Bare*^ 




D.-\-\j. 


01 


0.99 


99 


0.99 


1.0 


1.0 


1.0 


1.0 


1.0 


1.0 


1.0 
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09 


n 
u 


Dy 
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95 
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99 


n 


35 

O'J 


n 


35 

O'J 
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U.OO 


Al 





56 


56 





56 





19 





19 





19 





24 





24 





24 





32 





32 


n QO 

U.oZ 


A2 





55 


55 





55 





23 





23 





23 





28 





28 





28 





35 





36 


36 

\J •0\J 


A3 





52 


0.52 





52 





19 





19 





19 





24 





24 





24 





31 





31 


0.31 


A4 





49 


0.49 





49 





17 





17 
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21 





21 





21 





28 





28 


0.28 


A5 





52 


0.51 





52 





23 





23 





23 





27 





28 





28 





34 





34 


0.35 


A6 





48 


0.48 





48 





19 





19 





19 





24 





24 





24 





30 





30 


0.30 


A7 





46 


0.46 





45 





16 





16 





16 





20 





21 





21 





27 





27 


0.27 


A8 





50 


0.50 





50 





23 





23 





23 





28 





28 





28 





34 





34 


0.35 


A9 





47 


0.47 





47 





19 





19 





19 





24 





24 





24 





30 





30 


0.30 


AlO 





44 


0.44 





44 





16 





16 





16 





20 





20 





21 





26 





26 


0.27 


All 





47 


0.47 





47 





21 





21 





21 





25 





25 





25 





31 





31 


0.32 


A12 





48 


0.48 





48 





25 





25 





25 





29 





29 





29 





35 





35 


0.35 


Flat 





67 


0.67 





67 





68 





68 





68 





67 





67 





67 





62 





62 


0.62 



^Renormalized such that it is unity for the model with the highest marginal probability density distribution 
function value, for the data set and the appropriate beamwidth- and calibration-uncertainty treatment. 
When the likelihood normalization is set such that is it unity for Qrms-PS = 0, these highest marginal values 
are roughly 50, 1 • 10^, 6 • 10^ and 1 • lO^i for Ka2, Q2, Ka and Ka+Q respectively. 

'^Computed using eq. (8). Ignores beamwidth and calibration uncertainties. 

■^Accounts for beamwidth uncertainty using eq. (12). Ignores calibration uncertainty. 

■^Computed using eqs. (12) and (16). Accounts for beamwidth and calibration uncertainties. 
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Table 12: Numerical Values for Q,- 
Uncertainties^ 



_ps (in /iK), Accounting for Beamwidtli and Calibration 



Model (Oo, h, ^Bh^ 



Kal Ka2 Ka3 Ka4 



Kc 
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(0 
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75, 





0125) 


23 


1 


02 


(0 
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0175) 
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03 
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26 
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04 


(0 
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5 
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25 
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06 


(0 


3, 
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25 


8 
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(0 
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25 


8 


08 


(0 
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0175) 


23 


4 


09 


(0 
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23 


4 
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(0 
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23 


4 
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(0 


5, 
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5 
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(0 
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20 


5 
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(0 
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20 


5 
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(1 
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19 


8 


Al 


(0 
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6 


A2 


(0 
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23 
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A3 
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0175) 


23 


7 


A5 


(0 


3, 





70, 
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21 


3 


A9 


(0 
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21 


3 
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(0 
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0175) 


21 


3 
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(0 
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0125) 


20 


8 
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(1 
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20 


4 


Flat 














20 


4 



25.6 
20.5 
, 29.4 
' 23.6 
, 29.4 
' 23.6 
, 29.4 
' 23.6 
. 28.6 
' 23.0 
. 28.6 
' 23.0 
. 28.6 
' 23.0 
25.9 

- 20.8 
25.9 

- 20.8 
25.9 

- 20.8 
, 22.8 

18.3 

■ 22.8 
18.3 

, 22.8 
18.3 
. 22.0 
' 17.7 
• 28.4 
' 22.8 
'26.3 
21.1 
'26.3 
21.1 
'26.3 
21.1 
I 24.7 
' 19.8 
I 24.7 
' 19.8 
I 24.7 
' 19.8 
. 23.7 
' 19.0 
. 23.7 
' 19.0 
. 23.7 
' 19.0 
. 23.1 
' 18.5 
22.7 

■ 18.2 
22.6 

■ 18.2 



46 

45 
46 
47 
38 
39 
41 
32 
33 
34 
25 
26 
28 
23 
27 
26 
25 
23 
24 
23 
22 
24 
22 
21 
22 
23 
33 



66 
32 
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17 . 


1 45 

31 21 
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1 a 43 
16 ... 
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30 21 


65 
32 


17 . 
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67 

33 


1 ^ 45 

17 .: 


09 46 
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55 
26 


1 A 37 

14 ... 


38 

26 18 


57 
27 


1 /I 38 

14;. 


O'T 39 

27 ti 


59 
28 


1 c 39 

15 :. 


00 40 

28 19 


46 
22 


1 ri 30 

12 :. 


22^1 


47 
22 


1 O 31 
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00 32 

22 15 


49 
23 


1 o 33 

12 ... 
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37 
17 


r> o 24 

9.3:. 


1 25 
17^2 


38 
18 


n T 25 

9.7 .: 


1 26 
18^2 


40 
19 


1 n 27 

10:.. 


1 r> 27 
19l3 


33 
16 


O /I 22 

8.4... 


1 C 23 
I611 


39 
19 


n n 25 

9.9:.. 


1 n 27 
19 13 


38 
18 


9.7:.: 


1 26 
18^2 


35 
17 


9.1 2.^* 


17- 


34 
16 


8.722 


1623 


35 
17 


9.123 


172^ 


33 
16 


8.522 


16- 


31 
15 


8.021 


15?^ 


34 
16 


8.722 


16- 


32 
15 


8.221 


1 522 
i-J 10 


30 
14 


7.720 


1421 


32 
15 


8.221 


15?^ 


33 
16 


8.422 


16- 


48 
22 


12^2 


9932 

15 



2638 

ZU 

2638 

ZU 

9739 

^' 18 
9740 
^' 18 

2333 

23 ?3 
9434 

^^16 
10 27 
13 

19?^ 
20- 

15?^ 
16- 
16- 

14i°i 
16- 

16- 
1 521 
1420 

1 521 
1420 

138.8 
1420 

10 19 
^^"^8.9 
10 18 
^^"^8.4 
10 19 
1^*^8.9 
1420 
^^9.1 
27 
13 



14 

29 



19 



90 39 
^^20 

9-7 38 
^' 19 
90 39 
^^20 
90 40 
^^20 

2333 

24 
2535 

ZO 18 
10 27 

ly 14 
20- 

21^5 
16- 

16 23 
12 

1724 

^ ' 12 

1420 

10 

1723 
^ ' 12 

16- 

15 fi 
15 20 

1521 

1420 
10 
1419 

1 520 
10 

14 19 

10 18 
^^'^9.4 

1419 

1420 
10 
28 
14 



20 



^For each model, the first of the three entries in each of the last six columns is the value at which the 
probability density distribution function peaks. Ellipses as the lower entry in a vertical pair denotes a non- 
detection; the corresponding upper entry is the 2a (97.72% ET) upper limit. For detections, the vertical 
pair of numbers are the ±1(T (68.27% HPD) upper and lower limits. 

'^DMR values for open models is from two-year DMR (galactic-frame, quadrupole-excluded) normalization of 
Gorski et al. (1995) (also see RSBG), for fiat-A models from two-year DMR (ecliptic-frame, quadrupole- 
excluded) normalization of Bunn & Sugiyama (1995) (also see RSBG), and for fiat bandpower from two-year 
DMR (galactic-frame, quadrupole-excluded) normalization of Gorski et al. (1994). DMR la range accounts 
for DMR statistical and systematic uncertainties following Stompor et al. (1995) (also see RSBG). 
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Table 13: Numerical Values for Qrms-PS (in /^K), Accounting for Beamwidtli and Calibration 
Uncertainties^ 



Model 




(Oo, h, VtBh^) 


DMR 


Ql 


Q2 


Q3 


Q 


Ka + Q 


01 


(0 


1, 





75, 





0125) 


23 


1 25.6 
20.5 


25?^ 


"J" 26 


3828 


371° 


Of; 47 
'^■J 26 


02 


(0 


2, 





65, 





0175) 


26 


c 29.4 
23.6 


25 


of;48 
"J "J 26 


37-7 


36 f7 


3446 


03 


(0 


2, 





70, 





0125) 


26 


c 29.4 
"-^ 23.6 


25?^ 


"J" 26 


37^28 


37 1 


Of; 47 
'^■J 26 


04 


(0 
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75, 





0075) 


26 


c 29.4 

23 6 


26 f° 


0751 
1 27 
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3829 


36 f7 


05 


(0 


3, 





60, 





0175) 


25 


o 28.6 
23 


21 ?3 


30 i 


31 23 


31 23 


00 39 
■^^ 22 


06 


(0 


3, 





65, 





0125) 


25 


o 28.6 
O 23.0 


22^4 
12 


31 23 


09 44 
iiL 24 


32 f4 


30^3 


07 


(0 


3, 





70, 





0075) 


25 


o 28.6 
23.0 


99 35 
13 


09 44 
0^24 


33 


33 


01 42 
Oi 24 


08 


(0 


4, 





60, 





0175) 


23 


4 25.9 
^ 20.8 


17^8 


OK 35 
^0 19 




L\i ^g 


9432 
18 


09 


(0 


4, 





65, 





0125) 


23 


4 25.9 
20.8 


18?^ 


r,c36 


27 fo 


27 fo 


Of; 34 
^0 19 


OlO 


(0 


4, 





70, 





0075) 


23 


4 25.9 
20.8 




2737 
^' 20 


28^8 


28 f 


20 


Oil 


(0 


5, 





55, 





0175) 


20 


c 22.8 
18.3 


1422 


20^8 

15 


91 28 
15 


91 28 
16 


10 26 
15 


012 


(0 


5, 





60, 





0125) 


20 


c 22.8 
18.3 


1^8.2 


91 29 
16 


99 29 


99 29 


20 2^ 

15 


013 


(0 


5, 





65, 





0075) 


20 


c 22.8 
"-^ 18.3 


1 524 


9930 
16 


23 17 


23 ?? 


21 28 
^-L 16 


014 


(1 


0, 





50, 





0125) 


19 


o 22.0 
° 17.7 


10 19 
^•^ 7.0 


18?^ 


10 25 

ly 14 


19?^ 


1723 
^ ' 13 


Al 


(0 


1, 





90, 





0125) 


25 


a28A 
" 22.8 


^g36 


9930 


9930 


9930 


21 28 
^-L 16 


A2 


(0 


2, 





80, 





0075) 


23 


7 26.3 
' 21.1 


1422 


91 29 
16 


91 29 
16 


99 29 
16 


20 2^ 

15 


A3 


(0 


2, 





75, 





0125) 


23 


7 26.3 
' 21.1 


^533 


20 2^ 

15 


20 2^ 

15 


20 2^ 

15 


19 f5 


A4 


(0 


2, 





70, 





0175) 


23 


7 26.3 
' 21.1 


1431 


19?^ 


19 f4 


19?^ 


18^4 


A5 


(0 


3, 





70, 





0075) 


22 


9 24.7 
^ 19.8 


1021 
-L-J 7.4 


20 ?I 


20 ?I 


20 ?I 


19?^ 


A6 


(0 


3, 





65, 





0125) 


22 


9 24.7 
^ 19.8 


1431 


19?^ 


19?^ 


19?^ 


18^4 


A7 


(0 


3, 





60, 





0175) 


22 


9 24.7 
^ 19.8 


1329 


18?^ 


18^3 


18?! 


17?^ 


A8 


(0 


4, 





65, 





0075) 


21 


o 23.7 
19.0 


1020 
iO 7 2 


19?^ 


19 f4 


19?^ 


18?! 


A9 


(0 


4, 





60, 





0125) 


21 


o 23.7 
19.0 


ir,19 
-■-^6.7 


18?^ 


18^3 


18^4 


1723 

^ ' 13 


AlO 


(0 


4, 





55, 





0175) 


21 


o 23.7 
19.0 


^228 


17?^ 


17f3 


17f3 


1621 

12 


All 


(0 


5, 





60, 





0125) 


20 


o 23.1 
^18.5 


ir,19 
-■-^6.7 


18?^ 


18^3 


18^4 


1723 
^ ' 13 


A12 


(1 


0, 





50, 





0125) 


20 


422.7 
^ 18.2 


10 19 
^'^ 7.0 


18?^ 


1025 
1» 14 


1025 
1» 14 


1723 
^ ' 13 


Flat 














20 


422.6 
^ 18.2 


18?^ 


26?^ 


2737 
^' 20 


2737 
^' 20 


Of; 34 
^0 19 



^Conventions are the same as for Table 12. 
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Table 14: Numerical Values for 5Ti (in /iK), Accounting for Beamwidtli and Calibration 
Uncertainties^ 



Vlodel 






Kal 


Ka2 


Ka3 


Ka4 


Ql 


Q2 


q 


01 


(0 


1, 





75, 





0125) 


51 si 


1949 


35 23 


00 43 
20 


28^6 




42 


02 


(0 


2, 





65, 





0175) 


50 si 




34l°s 


on 43 
20 


28 fe 


40^^ 


42 


03 


(0 


2, 





70, 





0125) 


501 




34l°s 


00 43 
20 


28^6 


40 

^^30 


42 


04 


(0 


2, 





75, 





0075) 


501 


^g48 


34 fs 


00 43 
20 


28 fe 


40 

^^30 


42 


05 


(0 


3, 





60, 





0175) 


49 si 


18 5 


34 fs 


on 43 
20 


28 fe 


41 

30 


42 


06 


(0 


3, 





65, 





0125) 


491 


^g48 


34 fs 


on 43 
20 


28 fe 


40 

^^30 


42 


07 


(0 


3, 





70, 





0075) 


501 


^g48 


34l°s 


00 43 
20 


28 fe 


40 

^^30 


42 


08 


(0 


4, 





60, 





0175) 


481 


i8:'T 


34 fs 


on 43 
20 


28 fe 


41 

30 


42 


09 


(0 


4, 





65, 





0125) 


49 S3 


i8;'T 


34 fs 


00 43 
20 


28 fe 


41 

30 


42 


OlO 


(0 


4, 





70, 





0075) 


491 


^g48 


34 fs 


00 43 
20 


28 fe 


40 

^^30 


42 


Oil 


(0 


5, 





55, 





0175) 


481 


18 5 


34 f 


on 43 
20 


28 fe 


41 

30 


42 


012 


(0 


5, 





60, 





0125) 


481 


18 5 


34 f 


on 43 
20 


28 fe 


41 

30 


42 


013 


(0 


5, 





65, 





0075) 


49 S3 


i8:'T 


34 f 


00 43 
20 


28 f 


41 

30 


42 


014 


(1 


0, 





50, 





0125) 


48^^ 


^g46 


34 f 


30 fo 


28 f 


41 

30 


42 


Al 


(0 


1, 





90, 





0125) 


47^^ 


^g45 


34 f 


30 fo 


30^5 


40 58 


43 


A2 


(0 


2, 





80, 





0075) 


47^^ 


^g46 


34 f 


on 43 
20 


28 f 


4057 


42 


A3 


(0 


2, 





75, 





0125) 


47^^ 


1 46 

18 . 


34 f 


on 43 

30 


on 68 

30 _ 


42^1 


42 


A4 


(0 


2, 





70, 





0175) 


47^^ 


^g45 


34 f 


on 43 
20 


30 6« 


4958 


43 


A5 


(0 


3, 





70, 





0075) 


47^^ 


^g46 


34 f 


on 43 
20 


28 f 


4^57 


42 


A6 


(0 


3, 





65, 





0125) 


47^^ 


^g46 


34 f 


on 43 
20 


30^5 


4957 


42 


A7 
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^For each model, the first of the three entries in each of the last seven columns is the value at which the 
probability density distribution function peaks. Ellipses as the lower entry in a vertical pair denotes a non- 
detection; the corresponding upper entry is the 2a (97.72% ET) upper limit. For detections, the vertical 
pair of numbers are the ±1(T (68.27% HPD) upper and lower limits. 
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Table 15: Renormalized Maximum Values of the Probability Density Distribution Functions' 
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^Accounts for beamwidth and calibration uncertainties. Renormalized such that it is unity for the model 
with the highest maximum value of the probability density distribution function for the data set. The 
last line of the table gives this highest maximum likelihood value when the normalization is set such that 

i.((5rms-PS = IjK) = I. 
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Table 16: Renormalized Marginal Values of the Probability Density Distribution Functions' 
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^Accounts for beamwidth and calibration uncertainties. Renormalized such that it is unity for the model with 
the highest marginal probability density distribution function value for the data set. The last line of the 
table gives the marginal value for the model with highest marginal probability density distribution function 
value when the likelihoods are normalized such that ^(Qrms-PS = //K) = 1. 



FIGURE CAPTIONS 



Fig. 1. — CMB anisotropy multipole moments /(/ + l)C;/(27r) X 10^° (broken lines, scale on 
left function of multipole /, to / = 300, for selected models 01, Oil, 014, A2, AlO, 

and Flat, normalized to the two-year DMR maps. 01, Oil, 014, and Flat are normalized 
to the galactic-coordinates quadrupole-excluded maps, while A2 and AlO are normalized 
to the ecliptic-coordinates quadrupole-excluded ones. See Table 13 for model-parameter 
values. Also shown are two SP94 individual-channel, zero-lag, nominal beamwidth, window 
functions Wi (solid lines, scale on right axis), sensitive to the largest (Kal) and smallest 
(Q3) angular scales. See Table 1 for window function parameter values. 

Fig. 2. — Individual-channel, nominal beamwidth, zero-lag SP94 window functions Wi (solid 
lines), as a function of /, to / = 300. Also shown (hatched regions) is the uncertainty in 
the windows sensitive to the largest (Kal) and smallest (Q3) angular scales, due to the 
one standard deviation uncertainty in the beamwidths. See Table 1 for window function 
parameter values. The Wis are computed using eq. (1) of G95. The Wi used here assume 
that the SP94 beam is gaussian and do not explicitly account for the small ellipticity. 

Fig. 3. — {ST,-ms^)i (eq. [6]) as a function of /, to / = 250, for (upper panel) the larger 
beamwidth (upper Icr) Kal channel window, and for (lower panel) the smaller beamwidth 
(lower Icr) Q3 channel window, for the selected models shown in Fig. 1, normalized to the 
two-year DMR data. See Tables 2, 12, and 13 for numerical values. Note that the peak 
sensitivity of an SP94 individual-channel window function corresponds to a different angular 
scale in each of the models. 

Fig. 4. — SP94 data of G95. (a) Individual-channel Ka data, (b) Individual-channel Q data. 
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Fig. 5. — Likelihood functions as a function of Qrms-PS (to Qrms-PS = 60/iK) derived from 
analyses: ignoring beamwidtli and calibration uncertainties (solid lines); accounting for 
beamwidtli uncertainty but ignoring calibration uncertainty (dotted lines); and accounting 
for beamwidtli and calibration uncertainties (dashed lines). In (a) are those for the Ka2 
data (which does not have a detection), in (b) those for the Q2 data (which has one of the 
best individual-channel detections), in (c) those for the Ka data, and in (d) those for the 
Ka + Q data. The 6 individual panels shown for each of these SP94 data sets correspond 
to the 6 selected models of Fig. 1. See Table 2 for model-parameter values. In each panel 
the likelihood functions have been renormalized such that the peak value is unity for the 
one with the highest peak value. See Tables 8-16 for numerical values derived from the 
corresponding probability density distribution functions. 



Fig. 6. — Likelihood functions, as a function of Qrms-PS for the 6 selected models (01, Oil, 
014, A2, AlO, and Flat) of Fig. 1 (line styles are identical to those used in Figs. 1 and 
3), from the beamwidth- and calibration-uncertainty corrected analyses, for the Ka -|- Q, 
Ka, Ka2, and Q2 data sets. See Table 2 for model-parameter values, and Tables 8-16 for 
numerical values derived from the corresponding probability density distribution functions. 
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Fig. 7. — CMB anisotropy bandtemperature predictions and observational results. Hatched 
regions are what would be seen by a series of ideal, Kronecker-delta window-function, 
experiments, for the model normalized to the Icr range from the two-year DMR data (hatched 
region with smaller vertical extent), and for the model normalized to the Icr range from 
the SP94 combined Q (upper panel), Ka (central panel), and full Ka -|- Q (lower panel) 
data subsets (hatched region with larger vertical extent). These are evaluated using the 
appropriate model spectrum in eq. (7) for a Kronecker-delta Wi. Both normalizations 
account for all major known uncertainties. Also shown are the corresponding SP94 
individual-channel observational results. Solid squares correspond to detections, and are 
placed at the peak of the probability density distribution function. The vertical error 
bars on the points corresponding to detections are iblcr HPD error bars. Solid inverted 
triangles are 2cr ET upper limits corresponding to non-detections. From left to right, the 
observational data points correspond to the Kal to Q3 data results, and they are placed at 
/in determined from {ST,-ms^)i (eq. [6]). The horizontal lines on the SP94 individual-channel 
observational results represent the /-space width of the corresponding W;, and terminate at 
/e-0 5 determined from {ST,-ms^)i using the -\-lcr beamwidth Wi for the lower value of /g-0 5 
and the —Icr beamwidth Wi for the higher value. The 6 sets of 3 panels each [(a) to (f)] 
correspond to the 6 selected models of Fig. 1. See Table 2 for model-parameter values. 
Tables 12 and 13 for the SP94 combined data model normalizations and the two-year DMR 
model normalizations, and Table 14 for the SP94 individual-channel results. 
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